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Abstract. Here I review the current state of the field of optical stellar interferometry, 
concentrating on ground-based work although a brief report of space interferometry 
missions is included. We pause both to reflect on decades of immense progress in the 
field as well as to prepare for a new generation of large interferometers just now being 
commissioned (most notably, the CHARA, Keck and VLT Interferometers). First, 
this review summarizes the basic principles behind stellar interferometry needed by 
the lay-physicist and general astronomer to understand the scientific potential as well 
as technical challenges of interferometry. Next, the basic design principles of practical 
interferometers are discussed, using the experience of past and existing facilities to 
illustrate important points. Here there is significant discussion of current trends in the 
field, including the new facilities under construction and advanced technologies being 
debuted. This decade has seen the influence of stellar interferometry extend beyond 
classical regimes of stellar diameters and binary orbits to new areas such as mapping the 
accretion disks around young stars, novel calibration of the Cepheid Period-Luminosity 
relation, and imaging of stellar surfaces. The third section is devoted to the major 
scientific results from interferometry, grouped into natural categories reflecting these 
current developments. Lastly, I consider the future of interferometry, highlighting 
the kinds of new science promised by the interferometers coming on-line in the next 
few years. I also discuss the longer-term future of optical interferometry, including 
the prospects for space interferometry and the possibilities of large-scale ground-based 
projects. Critical technological developments are still needed to make these projects 
attractive and affordable. 
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1. Introduction 

This review will introduce the theory, technique, and scientific goals of optical 
stellar interferometry. By combining light collected by widely-separated telescopes, 
interferometrists can overcome the diffraction-limit of an individual telescope. The 
angular resolution achieved by current instruments is indeed astounding, < 5 x 10~^ 
radians (1 milli-arcsecond) , as are the engineering feats of maintaining sub-micron 
optical stability and coherence over hundreds of meters of pathlength while controlling 
polarization and dispersion over broad wavelength bandpasses. New capabilities are 
being applied in a wide variety of astrophysics contexts, including fundamental stellar 
parameters, novel ways to measure distances to stars, probing star formation and 
evolution, direct detection of extrasolar planets, and resolving cores of the nearest active 
galactic nuclei and brightest quasars. This current work is pointing the way towards 
next generation facilities, and I will close with a discussion of efforts to bring the angular 
resolution advantages of interferometers into space. 

1.1. Scope of Review 

The history of stellar interferometry spans more than a century, and a proper 
documentation of the rich history is beyond the scope of this review (see Lawson, 
2000a, for a historical overview of the field). The main purpose of this review will 

be to summarize the current state of the field of optical interferometry, including 
past scientific and engineering lessons, current astronomical motivations, and future 
goals and performance expectations. Interested readers may want to consult earlier 
reviews which tend to emphasize other topics, especially Shao and Colavita (1992a) 
and Quirrenbach (2001). In order to restrict the length, I will concentrate on long- 
baseline optical interferometry, giving only passing description to diffraction-limited 
single-aperture experiments (e.g.. Speckle Interferometry, Aperture Masking, Adaptive 
Optics). Further, I consider "optical interferometry" in a restricted sense to mean 
the light from the separate telescopes are brought together using optics, as opposed 
to heterodyne interferometry whereby the radiation at each telescope is coherently 
detected before interference (although I will discuss briefiy the important cases of the 
Intensity Interferometer and heterodyne interferometry using CO2 lasers). In practice 
then, "optical interferometry" is limited to visible and infrared wavelengths, and I will 
not discuss recent advances in mm-wave and sub-mm interferometry. 

1.2. The Organization of Review 

This review is divided into 4 major sections. The first reviews the basic theory behind 
optical interferometry and image reconstruction through a turbulent atmosphere. The 
second section explains the basic designs of interferometers and core modern technologies 
which make them work, including descriptions of current facilities. Major scientific 
results are outlined in the third section. The last section forecasts near-future science 
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potential as well as the long-term prospects of optical interferometry on the ground and 
in space. 

1.3. Nomenclature 

In this review, "optical" does not indicate the visible portion of the electromagnetic 
spectrum only, but generally refers to how the light is manipulated (using optics); in 
the context of interferometry, this will limit our discussion to wavelengths from the blue 
(~ 0.4/im) to the near-infrared (1-5/i) and mid-infrared (8-12yum). Typical angular 
units used in this paper are "milli-arcseconds" , or mas, where an arcsecond is the 
standard 1/3600 of a degree of angle. Astronomers often use the magnitude scale to 
discuss the wavelength-dependent flux density (power per unit area per unit bandwidth) 
from an astronomical source, where the bright star Vega {a Lyrae) is defined as mag 
(corresponding to a lOOOOK blackbody); the magnitude scale is logarithmic such that 
every factor of 10 brightness decrease corresponds to a flux "magnitude" increase of 2.5 
(e.g., a contrast ratio of 10 astronomical magnitudes is a factor of 10^). In addition, the 
unit Jansky (Jy) is often used for measuring the flux density of astronomical objects, 
1 Jy = 10-^^Wm-^Hz-\ 

A number of other basic astronomical units are used herein. The distance between 
the Earth and Sun is one astronomical unit, 1 AU ~ 1.5 x lO^^m. Stellar distances 
are given in units of parsecs (1 parsec is the distance to a star exhibiting a parallax 
angle of 1 arcsecond): written in terms of other common units of length, 1 pc 
~ 3.09 X 10^^m~ 3.26 light years. 

Lastly, I want to alert the reader (in advance) to Table 2 which will define all the 
interferometer acronyms used throughout the text. 

2. Basic Principles of Stellar Interferometry 

This section will review the basic principles of stellar interferometry. More detailed 
discussions of optical interferometry issues can be found in the recent published 
proceedings of the Michelson Summer School, "Principles of Long Baseline Stellar 

Interferometry" edited by P. Lawson (Lawson, 2000b), and the proceedings of the 2002 
Les Houches Eurowinter school edited by G. Pcrrin and F. Malbet (Pcrrin and Malbet, 
2002); earlier such collections also continue to play an important reference role (e.g., 
Perley ct al., 1986; Lagrange ct al., 1997). Other "classic" texts on the subjects of 
radio interferometry and optics are Thompson et al. (2001), Born and Wolf (1965), and 
Goodman (1985). 

2.1. Basics of Stellar Interferometry 

The basic principles behind stellar interferometry should be familiar to any physicist, 
founded on the wave properties of light as first observed by Thomas Young in 1803. 
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Figure 1. Young's two-slit interference experiment (monochromatic light) is presented 

to illustrate the basic principles behind stellar interferometry. On the left is the case 
for a single point source, while the case on the right is for a double source with the 
angular distance being half the fringe spacing. Note, the interference pattern shown 
represents the intensity distribution, not the electric field. 



This result is widely known through Young's "two-slit experiment," although two slits 
were not used in the original 1803 work. 

2.1.1. Young's Two-slit Experiment In the classical setup, monochromatic light from a 
distant ( "point" ) source impinges upon two slits, schematically shown in the left panel of 
Figure 1. The subsequent illumination pattern is projected onto a screen and a pattern 
of fringes is observed. This idealized model is realized in a practical interferometer 
by receiving light at two separate telescopes and bringing the light together in a beam- 
combination facility for interference (this will be discussed fully in §3; e.g., see Figure 10). 
The interference is, of course, due to the wave nature of light a la Huygens; the electric 
field at each slit (telescope) propagating to the screen with different relative path lengths, 
and hence alternately constructively and destructively interfering at different points 
along the screen. One can easily write down the condition for constructive interference; 
the fringe spatial frequency (fringes per unit angle) of the intensity distribution on the 
screen is proportional to the projected slit separation, or baseline 6, in units of the 
observing wavelength A (see Figure 1). That is. 

Fringe Spacing = A© = ^ radians (1) 



Fringe Spatial Frequency = u—^ radians"^ (2) 

Imagine another point source of light (of equal brightness, but incoherent with 
the first) located at an angle of A/ (26) from the first source (see right of panel of 
Figure 1). The two illumination patterns are out of phase with one another by 180°, 
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hence cancelling each other out and presenting a uniformly illuminated screen. Clearly 
such an interfering device (an "interferometer" ) can be useful in studying the brightness 
distribution of a distant "stellar" object. This application of interferometry was first 
proposed by Fizeau in 1868 (Fizeau, 1868) and successfully applied by Michclson to 
measure the angular diameters of Jupiter's moons (Michelson, 1890, 1891) in 1891 and 
later (with Pease in 1921) to measure the first angular size of a star beyond the Sun 
(Michelson and Pease, 1921) (see §3.1 for further details on the early history of optical 
interferometry) . 

2.1.2. Angular Resolution The ability to discern the two components of a binary star 
system is often used to gauge the spatial resolution of an instrument, be it a conventional 
imaging telescope or a separated-element interferometer. Classical diffraction theory has 
established the "Rayleigh Criterion" for defining the (diffraction-limited) resolution of 
a filled circular aperture of diameter D: 

Resolution of Telescope = AOteiescope = 1.22 -^radians (3) 

This criterion corresponds to the angular separation on the sky when one stellar 
component is centered on the first null in the diffraction pattern of the other; the binary 
is then said to be resolved. A similar criterion can be defined for an interferometer: an 
equal brightness binary is resolved by an interferometer if the fringe contrast goes to 
zero at the longest baseline. As motivated in the last paragraph, this occurs when the 
angular separation is ^, where b is the baseline. Hence, 

Resolution of Interferometer = AOinterferometer = --rradians (4) 

While these two criteria are somewhat arbitrary, they are useful for estimating the 
angular resolution of an optical system and are in widespread use by the astronomical 
community. 



2.1.3. Complex Visibility One can be more quantitative in interpreting the fringe 
patterns observed with an interferometer. The fringe contrast is historically called the 
visibility and, for the simple (two-slit) interferometer considered here, can be written as 
y _ /max - /min _ Fringe Amplitude . , 

/max + /min Average Intensity 
where /max and /min denote the maximum and minimum intensity of the fringes. Hence, 
the left and right fringe patterns of Figure 1 have visibilities of one and zero respectively. 

The Van Cittcrt-Zcrnike Theorem (see Thompson et al., 2001, for complete 
discussion and proof) relates the contrast of an interferometer's fringes to a unique 
Fourier component of the impinging brightness distribution. In fact, the visibility is 
exactly proportional to the amplitude of the image Fourier component corresponding to 
the (spatial) fringe spatial frequency defined above {u — b/X rad~^). Also, the phase of 
the fringe pattern is equal to the Fourier phase of the same spatial frequency component. 
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Figure 2. This figure shows simple one-dimensional images and their corresponding 
visibility curves. The left panels are the images while the right panels correspond to 

the Fourier amplitudes, i.e. the visibility amplitudes. Note that "large" structure in 
image-space result in "small" structure in visibility-space. 



The Van Cittert-Zernike Theorem can be expressed concisely in mathematical 
terms. Consider that the astronomical target emits light at frequency v over only a 
very small portion of the sky with specific intensity Ii,{9, (p), so small that the spherical 
coordinates Oq + S6 and (f)Q + 6(f) can be interpreted as Cartesian coordinates x^i and ijq 
centered around 0o on the plane of the sky. We can write the interferometer 

response (amplitude and phase of the fringes) as the frequency-dependent complex 
visibility Vi,{u, v), defined as the Fourier Transform of the brightness distribution It,{fa), 
normalized so that — 0) — 1. 



„ / dxn dyn 4(ra)e ^'''( ^ '^") 
IV.lyjle-^^-" = '-^ (6) 



sn 



Total Specific Flux 

using the following notation: 

rh = {xn,yn) 
D 

— = the basehne vector D projected onto the plane 
of the sky in units of wavelength A 
= (m, v) [Common Notation] (7) 

Figure 2 shows some simple examples of one-dimensional images and the 
corresponding visibility curves. The top panels show the case of an equal binary 
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system, where both components are unresolved. The periodicity in the visibihty-space 
corresponds to the binary separation. The middle set of panels is representative of a 
compact, but resolved, source (such as a star surrounded by an optically-thick dust 
shell). The small image size means there is more high spatial frequency information, 
and this is why the corresponding visibility curve is non-zero even at high resolution. 
Lastly, the bottom panels show an image of a unresolved star (with 10% of the total flux) 
surrounded by larger-scale structure (this is expected when a star is surrounded by an 
optically-thin envelope of dust) . The large-scale structure (containing 90% of the total 
flux) can be seen to be "resolved" on short basehnes (at low spatial frequency), while 
the point-source remains unresolved out to the highest spatial frequency. Note that the 
visibility plateaus at 0.10, corresponding to the fraction of the total flux which is left 
unresolved. This is easy to understand since the Fourier Transform is linear; that is, the 
(complex) visibility of a point-source and extended structure is equal to the visibility of 
the point-source plus the visibility of the extended structure separately. This property 
of linearity is very helpful in interpreting simple visibility curves. 

Most astronomical objects are not one-dimensional, and the two-dimensional space 
of spatial frequencies is called the Fourier Plane, or the (u,v) plane, named after the (u,v) 
coordinates deflned in Eq.7. Further, in general we must consider both the visibility 
amplitude and the visibility phase. For example, consider the equal binary system 
depicted in Figure 3. The complex visibility can be easily written by choosing the origin 
midway between the two components. Note the abrupt phase jump when the visibility 
amplitude goes through a null. These discontinuities are smoothed out when the two 
components are not precisely equal. 

2.2. Atmospheric Problems 

An incoming plane wave from a stellar source is corrupted as it propagates through the 
turbulent atmosphere. Variations in the column density of air along different paths cause 
the effective pathlength to vary, introducing wavefront distortion. If these distortions 
become a signiflcant fraction of a wavelength across the aperture of a telescope, the size 
of the image formed will not be diffraction-limited by the primary mirror, but rather by 
the coherence scale of the incoming wavefront. The transverse distance over which one 
expects rms pathlength difference to be A/2.4 has been deflned as the Fried parameter 
and is denoted by ro(A) (Fried, 1965); hence telescope apertures larger than ro(A) can 
expect significant degradation of image quality (when observing at wavelength A) due 
to atmospheric effects. In fact, for an ro diameter circular patch, the rms phase error is 
~1.03 radians. At A =500 nm, ro is typically 10 cm (toward zenith) at average observing 
sites and hence even small telescopes can not be used at their diffraction limit in the 
visible. In such cases, the observed angular size of a point source will be determined 
entirely by ro(A) at a given wavelength, and is known as the seeing disk size, ©seeing (A). 
The Kolmogorov theory of turbulence (Kolmogorov, 1961) predicts that ro(A) oc \^/^, 
and hence the seeing size, ©seeing (A) oc oc A is only weakly dependent on the 



Optical Interferometry in Astronomy 

Example:An Equal Binary 



phase 
center 



separation p 



Visibility Amplitude 




-40 -20 20 40 
u 



Visibility Phase 



40 
20 


-20 
-40 




-40 -20 20 40 
u 



Figure 3. This figure shows the complex visibility for an equal binary system in 

the 2-dimension (u,v) plane. With the above choice for the phase center, the Fourier 
phases can be represented simply. Notice the abrupt phase jumps when the visibility 
amplitude goes through a null. This figure is reproduced through the courtesy of 
the NASA/ Jet Propulsion Laboratory, California Institute of Technology, Pasadena, 
Cahfornia (Monnier, 2000). 



wavelength (Fried, 1965). An example of the phase delays associated with a snapshot of 
Kolmogorov turbulence can be seen in Figure 4 for 12-m square, corresponding roughly 
to the size of the largest telescopes today (e.g., the Keck telescopes). 

Another consequence of turbulence is that the image distortion varies across the 
sky, although stars located close together suffer similar seeing effects. The angle over 
which image distortions are correlated is called the "isoplanatic" angle, and is only 
a few arcseconds in the visible and about an arcminute in the near-infrared. This 
angle is determined by the vertical distribution of the turbulence - obviously low-level 
turbulence would induce correlated image distortions over larger sky angles than the 
same turbulent layer located higher up in the atmosphere. The isoplanatic angle is a 
critical parameter for the field-of-view of adaptive optics systems which actively sense 
and correct for atmospheric turbulence in realtime using a deformable mirror. 

Another important atmospheric diagnostic is the coherence time, to- Typically, one 
assumes a "frozen" turbulence model in which the atmospheric density perturbations 
are assumed constant over the time it takes wind to blow them across a given aperture 
(also known as Taylor's hypothesis of frozen turbulence). This motivates a convenient 
estimate for the coherence time: to(^) = '''o(-^)/^wmd) where fwmd is the wind speed. At 
most sites, wind speeds are ~10ms^^ and so to ~ 10 ms at 500 nm. Further discussion 
of atmospheric turbulence and degradation of astronomical images can be found in 
Kolmogorov (1961), Roddier (1981), and Roddier et al. (1982). 
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Figure 4. This figure shows a typical reahzation of Kolmogorov turbulence (ro = 

50 cm at A = 2.2 /im); each solid contour line represents ^ of wavefront distortion. 
Some areas of the aperture show coherent areas larger than ro, and some much smaller; 
ro is a statistical property of atmospheric turbulence and wavefront perturbations occur 
over a wide range of scales. 



These parameters, ro(A) and to(A), are extremely important for the design of an 
interferometer, because the value of ro hmits the useful size of the collecting aperture, 
while to limits the coherent integration time. Both of these are crucial for predicting an 
interferometer sensitivity to faint objects and much debate surrounds the best estimates 
for these parameters at various sites (e.g., Dyck and Howell, 1983; Roddier et al., 1990; 
ten Brummelaar, 1994; Treuhaft et al., 1995). This topic is revisited in §3.4 when I 
discuss the limiting magnitude of current interferometers. 

It is well-known that ro(A) and to (A) depend greatly on the observing site, and we 
now consider the unique seeing conditions of Mauna Kca, Hawaii, as an example. The 
coherence scale is unusually long due to the highly laminar flow of the Pacific winds over 
the peak of the mountain (elevation 4200 m); ro usually lies between 10 and 40 cm at 
500 nm (Wizinowich, 1999). However, the fast winds of the overhead jet stream result 
in very short coherence times: to between 1.5 and 10 ms (Wizinowich, 1999). It should 
be emphasized that seeing is notoriously difficult to characterize due to large variations 
in time (both on short time scales as well as seasonal ones) as evidenced by the large 
range of ro(A = 500 nm) and to(A = 500 nm) values just given. 

2.2.1. Atmospheric Phase Errors The fluctuating amount of integrated atmospheric 
pathlength above each telescope introduce wavefront time delays which show up as phase 
shifts in the measured fringes in an interferometer, as illustrated in Figure 5. In this 
figure, an optical interferometer is represented again by a Young's two-slit experiment, as 
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Figure 5. Atmospheric time delays or phase errors at telescopes cause fringe shifts, 
as can be seen through analogy with Young's two-slit experiment. This figure is 
reproduced through the courtesy of the NASA/ Jet Propulsion Laboratory, California 
Institute of Technology, Pasadena, California (Monnier, 2000). 



discussed earlier in this section. The spatial frequency of these fringes is determined by 
the distance between the slits (in units of the wavelength of the illuminating radiation). 
However if the pathlength above one slit is changed (due to a pocket of warm air moving 
across the aperture, for example) , the interference pattern will be shifted by an amount 

depending on the difference in pathlength of the two legs in this simple interferometer. 
If the extra pathlength is half the wavelength, the fringe pattern will shift by half a 
fringe, or vr radians. The phase shift is completely independent of the slit (telescope) 
separation, and only depends on slit-specific (telescope-specific) phase delays. 

The most obvious impact of atmospheric phase delays is that the assumptions of the 
van Cittert-Zernike theorem no longer apply, and that the measure fringe phase can no 
longer be associated with the Fourier phase of the sky brightness distribution (the fringe 
amplitude retains its original meaning, since phase changes do not change the measured 
fringe amplitudes for short exposures). The corruption of this phase information has 
serious consequences, since imaging of non-centrosymmetric objects rely on the Fourier 
phase information encoded in this intrinsic phase of interferometer fringes. Without this 
information, imaging can not be done except for simple objects such as disks or round 
stars. Fortunately, a number of strategies have evolved to circumvent these difficulties. 

2.2.2. Phase Referencing Possible methods for recovering this phase information using 
phase referencing techniques are discussed in Chapter 9 (written by A. Quirrenbach) 
of "Principles of Long Basehne Stellar Interferometry (Lawson, 2000b). Few scientific 
results have resulted from phase referencing techniques to date, but this is expected to 
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change over the coming decade as sophisticated new instruments are being deployed. 
Here, I mention a few of the most promising methods: 

(i) Neeirby sources. If a bright point reference source (or source with well-known 
structure) lies within an isoplanatic patch (see Quirrenbach, 2000), then its fringes 
will act as a probe of the atmospheric conditions. By measuring the instantaneous 
phases of fringes from the bright reference source, one can correct the corrupted 
phases on a ncaby "target" source. This has been applied to narrow-angle 
astromctry where fringe phase information is used for determining precise relative 
positions of nearby stars (Shao and Colavita, 1992b; Colavita et al., 1999; Lane 
et al., 2000a); see Figure 31 for some preliminary results published by the Palomar 
Testbed Interferometer. While it would be very valuable to use an artificial guide 
star for phase-referencing a long baseline interferometer, current laser beacons are 
too spatially extended (resolved) to produce interferometric fringes. 

(ii) A$ Monitoring. In the miUimeter and sub-millimeter, phase shifts caused by 
fluctuations in atmospheric water vapor column density can be monitored by 
observing its line emission. This information can be used to phase-compensate 
the interferometer, allowing longer coherent integrations and accurate fringe phase 
determination on the target (Wiedner 1998, and references therein). In the mid- 
infrared, strategies to actively monitor ground-level turbulence using temperature 
sensors are being explored by the Infrared Spatial Interferometer group (Short et al., 
2002) at Mt. Wilson motivated by recent atmospheric studies (e.g., Bester et al. 
1992). Townes (2002) recently proposed that realtime- monitoring of Rayleigh or 
Raman backscattering might be used to correct for atmospheric column density 
variations in the context of optical interferometers, but this method has not yet 
been validated. 

(iii) Multi-wavelength. Another possibility is to observe a target at multiple 
wavelengths and to use data from one part of the spectrum to calibrate another. 
For example, one might use fringes formed by the continuum emission to phase 
reference a spectral line (e.g., Vakili et al. 1997). To use this method, one must 
assume knowledge about the brightness distribution at one of the wavelengths being 
used. 

Currently, phase referencing is not possible with most current beam combiners in 
operation, either due to low spectral resolution or limited field-of-view. In order to 
recover phase information, one must one must make use of the closure phases. 

2.2.3. Closure Phases Consider Figure 6 in which a phase delay is introduced above 
Telescope 2 of a 3-telescopc array. As discussed in the last section, this additional 
delay causes a phase shift in the fringe detected between telescopes 1-2. Note that a 
phase shift is also induced for fringes between telescopes 2-3; however, this phase shift is 
equal and opposite to the one for telescopes 1-2. Hence, the sum of three fringe phases, 
between 1-2, 2-3, and 3-1, is insensitive to the phase delay above telescope 2. This 
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Figure 6. This figure explains the principle behind closure phase analysis. Phase 
errors introduced at any telescope causes equal but opposite phase shifts, canceling 
out in the closure phase (figure after Rcadhcad et al. 1988). This figure is reproduced 
through the courtesy of the NASA/ Jet Propulsion Laboratory, California Institute of 
Technology, Pasadena, California (Monnier, 2000). 



argument holds for arbitrary phase delays above any of the three telescopes. In general, 
the sum of three phases around a closed triangle of baselines, the closure phase, is a 
good interfcromctric observable; that is, it is independent of telescope-specific phase 
shifts induced by the atmosphere or optics. 

The idea of closure phase was first introduced to compensate for poor phase 
stability in early radio VLBl work (Jennison 1958). Application at higher frequencies 
was first mentioned by Rogstad (1968), but only much later carried out in the 
visible/infrared through aperture masking experiments (Baldwin et al. 1986a; Haniff 
et al. 1987a; Readhead et al. 1988; Haniff et al. 1989). Currently three separate- 
element interferometers have succeeded in obtaining closure phase measurements in the 
visible/infrared, first at COAST (Baldwin ct al. 1996a), soon after at NPOI (Benson 
et al. 1997a), and most recently at IOTA (Traub, 2002). 

How can these closure phases be used to figure out the Fourier phases which are 
needed to allow an image to be reconstructed? Each closure triangle phase can be 
thought of as a single linear equation relating three different Fourier phases (assuming 
none of the baselines are identical), which we desire to solve for; hence, we must count 
the number of linear equations available and compare to the number of unknowns. For 
N telescopes, there are "N choose 3," (^^^ — ^ possible closing triangles. 

However, there are only ( ^) = (^)(^~^) independent Fourier phases; clearly not all the 



,2/ 2 

closure phases can be independent. The number of independent closure phases is only 
(^2^^) ~ i^tzMltz^^ equivalent to holding one telescope fixed and forming all possible 
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Table 1. Phase information contained in the closure phases alone 



Number of Number of Number of 

Telescopes Fourier Phases Closing Triangles 



Number of Independent 
Closure Phases 



Percentage of 
Phase Information 



3 3 1 

7 21 35 

21 210 1330 

27 351 2925 

50 1225 19600 



1 

15 
190 

325 
1176 



33% 
71% 
90% 
93% 
96% 



triangles with that telescope (as discussed by Readhead et al., 1988). The number of 
independent closure phases is always less than the number of phases one would like to 
determine, but the percent of phase information retained by the closure phases improves 
as the number of telescopes in the array increases. Table 1 lists the number of Fourier 
phases, closing triangles, independent closure phases, and recovered percentage of phase 
information for telescope arrays of 3 to 50 elements. For example, approximately 90% 
of the phase information is recovered with a 21 telescope interferometric array (e.g., 
Readhead et al. 1988). As discussed in the next section, this phase information can be 
coupled with other image constraints (e.g., finite size and positivity) to reconstruct the 
source brightness distribution. 

In addition to the mathematical (linear algebra) interpretation of closure phases, 
there are a few other important properties worth noting. 

• For sources with point- symmetry (otherwise known as centro- symmetry), all the 
closure phases are either 0° or 180°. It is easy to prove this by imagining the 
image-center ( "phase-center" ) at the point of centro-symmetry. 

• Closure phases are not sensitive to an overall translation of image. A translation is 
indistinguishable from atmospheric phase delays for any given closing triangle. 

• The closure phases are independent of telescope-specific phase errors, however non- 
zero closure phases from a point source can result from having non-closing triangles 
and phase delays in the beam combiner (e.g., for a three-telescope pair- wise beam 
combiner) . 

2.3. Image Reconstruction 

While very few images have been made by today's optical interferometers, new telescope 
arrays are now being commissioned which will make true imaging interferometry 
straightforward. Because these new imaging capabilities are likely to have significant 
impacts over the next decade, I wish to review the basic principles of apertures synthesis 
imaging. However, I will restrain myself from excessive elaboration here, and instead 
refer the interested reader to the extensive radio interferometry literature, especially 
regarding "Very Long Baseline Interferometry (VLBI)." 
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While modeling visibility and closure phase data with simple models is useful, one 
would like to make an image unbiased by theoretical expectations. Since any image can 
be alternatively represented by its Fourier components, the collection of all "interesting" 
components can allow the interferometric data to be inverted, thus reconstructing an 
estimation of the image brightness distribution. The collection of a large number 
of Fourier components is greatly aided by increasing the number of telescopes, since 
independent combinations of telescopes increase with the number of telescopes to the 
second power ■ gge last section) 

With a large number of measurements, images of arbitrary complexity should 
be attainable using visible/infrared interferometers and reliable closure phase 
measurements. The importance of "filling up" the (u,v) plane with measurements when 
imaging is discussed more fully in §2.4.3. The next subsections will discuss strategies 
currently employed based on the techniques of VLBI in the radio. 

2.3.1. Guiding Principles The goals of an image reconstruction procedure can be 
stated quite simply: find an image which fits both the visibility amplitudes and closure 
phases within experimental uncertainties. However in practice, there are an infinite 
number of candidate images which satisfy these criteria, because interferometric data 
is always incomplete and noisy. Furthermore, the closure phases can not be used to 
unambiguously arrive at Fourier phase estimates as stated above, even under ideal noise- 
free conditions. 

Additional constraints must be imposed to "select" an image as the best-estimate 
of the true brightness distribution (to "regularize" this ill-posed inverse problem). These 
constraints introduce correlations in the Fourier amplitudes and phases, and essentially 
remove degrees of freedom from our inversion problem. Some of the most common (and 
reasonable) constraints are described below. 

• Limited Field-of-View. This constraint is always imposed in aperture synthesis 
imaging, even for a fully-phased array (e.g., VLA). Limiting the field-of-view 
introduces correlations in the complex visibility in the (u,v) plane. This is 
a consequence of the Convolution Theorem, a multiplication in image-space is 
equivalent to a convolution in the corresponding Fourier-space. 

• Positive-Definite. Since brightness distributions can not be negative, this 
is a sensible constraint (although not appropriate in some cases, such as 
for reconstruction of Stokes/polarization components or imaging spectral line 
absorption). While clearly limiting the range of "allowed" complex visibilities, 
there are few obvious, intuitive effects in the Fourier-plane; one is that the visibility 
amplitude is maximum at zero spatial frequency. The Maximum Entropy Method 
(see §2.3.3) naturally incorporates this constraint. 

• "Smoothness." The Maximum Entropy Method (MEM), for instance, selects the 
"smoothest" image consistent with the data. See §2.3.3 for more discussion of 
MEM. 
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• a priori Information. One can incorporate previously known information to 
constrain the possible image reconstructions. For instance, a low resolution image 
may be available from a single-dish telescope. Another commonly encountered 
example is point source embedded in nebulosity; one might want the reconstruction 
algorithm to take into account that the source at the center is point-like from 
theoretical arguments. 

For a phased interferometric array (e.g., the Very Large Array) where the Fourier 
phases are directly measured (avoiding the need for closure phases), one can use a 
number of aperture synthesis techniques to produce an estimate of an image based 
on sparsely sampled Fourier components. These procedures basically remove artifacts, 
i.e. sidelobes, of the interferometer's point-source response arising from uncomplete 
sampling of the (u,v) plane. These procedures do not incorporate closure phases, 
but work by inverting the Fourier amplitudes and phases to make an image. A 
brief explanation of the popular algorithms CLEAN and MEM follow with additional 
references for the interested reader. See Perley et al. (1986) for essays on these topics 
aimed at radio astronomers. 

2. 3. 2. CLEAN Originally described by Hogbom ( 1974) , CLEAN has been traditionally 
the most popular algorithm for image reconstruction in the radio because it is both 
computationally efficient and intuitively understandable. Given a set of visibility 
amplitiidcs and phases over a finite region of the Fourier plane, the "true" image can 
be estimated by simply setting all other spatial frequencies to zero and taking the 
(inverse) Fourier Transform. As one might expect, this process leads to a whole host 
of image artifacts, most damaging being positive and negative "sidelobes" resulting 
from non-complete coverage of the Fourier plane; we call this the "dirty map." The 
unevenly-filled Fourier plane can be thought of as a product of a completely-sampled 
Fourier plane (which we desire to determine) and a spatial frequency mask which is 
equal to 1 where we have data and elsewhere. Since multiplication in Fourier space is 
identical to convolution in image space, we can take the Fourier transform of the spatial 
frequency mask to find this convolving function; we call this the "dirty beam." Now the 
image reconstruction problem can be recast as a "deconvolution" of the dirty map with 
the dirty beam. 

The dirty map is CLEANed by subtracting the dirty beam (scaled to some fraction 
of the map peak) from the brightest spot in the dirty map. This removes sidelobe 
structure and artifacts from the dirty map. Repeating this process with dirty beams of 
ever decreasing amplitudes leads to a series of delta-functions which, when combined, fit 
the interferometric data. For visualization, this map of point sources is convolved with 
a Gaussian function whose FWHM values are the same as the dirty beam; this removes 
high spatial resolution information beyond the classic "Rayleigh" criterion cutoff. One 
major weakness with CLEAN is that this smoothing changes the visibility amplitudes, 
hence the CLEANed image no longer strictly fits the interferometric data, especially 
the spatial frequency information near the diffraction limit. Another weakness is that 



Optical Interferometry in Astronomy 



16 



CLEAN does not directly use the known uncertainties in the visibihty data, and hence 
there is no natural method to weight the high SNR data more than the low SNR data 
during image reconstruction. Further discussion of various implementations of CLEAN 
can be found in Clark (1980), Schwab (1984), Cornwell (1983), and Chapter 7 of Perley 
et al. (1986) by T. Cornwell. 

2.3.3. MEM The maximum entropy method (MEM) makes better use of the highest 
spatial frequency information by finding the smoothest image consistent with the 
interferometric data. While enforcing positivity and conserving the total flux in the 
frame, "smoothness" is estimated here by a global scalar quantity S. the "entropy." If /j 
is the fraction of the total flux in pixel i, then S — — Z^j /i In y- after the thermodynamic 
quantity; is known as the image prior and must be specified by the user. The MEM 
map fi will tend toward Jj when there is little (or noisy) data to constrain the fit. Often 
li is assumed to be a uniformly bright background, however one can use other image 
priors if additional information is available, such as the overall size of the source which 
may be known from previous observations. 

Mathematically, MEM solves the multi-dimensional (N=number of pixels) 
constrained minimization problem which only recently has become computationally 
realizable on desktop computers. Maintaining an adequate fit to the data (Sx^ ~ 
number of degrees of freedom), MEM reconstructs an image with maximum S. MEM 
image reconstructions always contain some spatial frequency information beyond the 
diffraction limit in order to keep the image as "smooth" as possible consistent with the 
data. Because of this, images typically have maximum spatial resolution a few times 
smaller than the typical Rayieigh-type resolution encountered with CLEAN ("super- 
resolution" ) . Further discussions of MEM and related Bayesian methods can be found 
in Pina and Puetter (1992), Narayan and Nityananda (1986), SkiUing and Bryan (1984), 
Gull and SkiUing (1983), and Sivia (1987). 

Unfortunately, MEM images also suffer from some characteristic artifacts and 
biases. Photometry of MEM-deconvolved images is necessarily biased because of the 
positivity constraint; any noise or uncertainty in the imaging appears in the background 
of the reconstruction instead of the source, systematically lowering the estimated fluxes 
of compact sources. Also, fields containing a point source embedded in extended 
emission often show structure reminiscent of Airy rings, the location of the rings being 
influenced by the wavelength of the observation and not inherent to the astrophysical 
source. Fortunately, these imaging artifacts are greatly alleviated for asymmetric 
structures, when closure phases and not the visibility amplitudes play a dominant role 
in shaping the reconstructed morphology. 

2.3.4- Including Closure Phase Information The above algorithms were originally 
designed to use Fourier amplitudes and phases, not closure phases. In order to use these 
algorithms, one has to come up a way to estimate the Fourier phases, when only the 
closure phases are available. Early image reconstruction algorithms incorporated closure 
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Figure 7. This is a flow diagram for a incorporating closure phase information into 
CLEAN/MEM aperture synthesis imaging algorithms based on the "self-calibration" 
procedure of Cornwell and Wilkinson (1981). This figure is reproduced through the 
courtesy of the NASA/ Jet Propulsion Laboratory, California Institute of Technology, 
Pasadena, California (Monnier, 2000). 



phase information by using an iterative scheme (Thompson et al. 1986; Readhead and 
Wilkinson 1978). The following steps summarize this process: 

(i) Start with a Fourier "phase model" based on either prior information or setting all 
phases to zero. 

(ii) Determine candidate phases by using some values from the "phase model" and 
enforcing all the (self-consistent) closure phase relations (see §2.2.3). 

(iii) Using CLEAN or MEM, perform aperture synthesis mapping on the given 
visibilities and candidate phases. At this stage, image constraints such as positivity 
and/or finite support are applied. 

(iv) Use this image as a basis for a new "phase model." 

(v) Go to step 2 and repeat until the process converges to a stable image solution. 

Cornwell and Wilkinson (1981) introduced a modification of the above scheme by 
explicitly solving for the telescope-specific errors as part of the reconstruction step. 
Hence the measured (corrupted) Fourier phases are fit using a combination of intrinsic 
phases (which are used for imaging using CLEAN/MEM) plus telescope phase errors. 
In this scheme, the closure phases arc not explicitly fit, but rather are conserved in 
the procedure since varying telescope-specific errors can not change any of the closure 
phases. Figure 7 shows a fiow diagram for this procedure; thoughtful consideration 
is required in order to fully understand the power and elegance of self-calibration, 
affectionately known as "self-cal." 
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Self-calibration works remarkably well for large number of telescopes, but requires 
reasonably high signal-to-noise ratio (SNR>5) complex visibilities. Once the SNR 
decreases below this point, the method completely fails. The conceptualization of 
solving for telescope-specific errors, while useful for the radio, is not applicable for 
visible/infrared interferometry where the good observables are the closure phases 
themselves, not corrupted Fourier phases. This is because the time-scale for phase 
variations in the visible/infrared is much less than a second, as opposed to minutes/hours 
in the radio. 

Of course, the self-cahbration iteration loop can be sidestepped altogether by fitting 
directly to all the data, the visibility amplitudes and closure phases, using MEM or some 
other regularization scheme. This would have the added advantage of allowing all the 
measurement errors to be properly addressed, theoretically resulting in the optimal 
image reconstruction. Buscher (1994) suggested this approach, but there has been little 
demonstrated progress in this method to date. 1 anticipate revived activity as more 
interferometers with "imaging" capability begin to produce data. 

2.3.5. Speckle Interferometry Another related interferometric technique which permits 
diffraction-limited observation through a turbulent atmosphere using a single fiUed- 

apcrturc telescope is "speckle interferometry," the promise of which was first realized by 
Labeyrie (1970a) in 1970. In §2.2, 1 claimed that observed angular size of a point source 
will be determined entirely by ro at a given wavelength, and is known as the seeing 
disk size, Osccing- However, this is only true for a long-exposure image. A single short- 
exposure image of a star actually consists of a network of small "speckles" extending 

over O seeing* 

In the original formulation of speckle interferometry, short exposures of an 
astrophysical object are made to freeze this "speckling" induced by the turbulent 
atmosphere. The amount of high-resolution structure in the speckle pattern, as 
quantified by its power spectrum, is a measure of two things: 1) the quality of the 
atmospheric seeing, and 2) the high resolution structure in the object of interest. 
Observing a nearby point-source star allows the calibration of the seeing contribution 
and thus the extraction of interferometric visibility measurements out to the diffraction 
limit of the telescope (i.e., the longest baseline). In analyzing this situation, one can 
think of many virtual subapertures (with size equal to the coherence length ro) spread 
across the full telescope, with fringes forming between all the subaperture pairs. After 
the original formulation by Labeyrie, it was discovered that the Fourier phases could 
also be estimated from such data (e.g., Knox and Thompson, 1974; Weigelt, 1977). 

Speckle interferometry data is often reduced using the "bispectrum," which permits 
a direct inversion from the estimated Fourier amplitudes and phases. The bispectrum 
Bijk = VijVjkVki is formed through triple products of the complex visibilities around 
a closed triangle, where ijk specifies the three aperture locations on the pupil of the 
telescope. One can see the bispectrum is a complex quantity, and that the bispectrum 
phase is identical to the closure phase. Interestingly, the use of the bispectrum for 
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reconstructing diffraction-limited images was developed independently (Wcigelt 1977; 
Hofmann and Weigelt 1993) of the closure phase techniques, and the connection between 
the approaches realized only later (Roddier, 1986; Cornwell, 1987). 



2.4- Other Important Considerations 

2.4- 1- Coherency The tolerance for matching pathlengths in an interferometer depend 
on the desired spectral bandwidth. In the limit of monochromatic light, such as for 
a laser, interference will occur even when pathlengths of an interferometer are highly 
mismatched. For broadband ("white") light, the number of fringes in an interferogram 
is equal to the inverse of the fractional bandwidth: A^fringes ~ -^x- Hence, for broad 
band observations (~20% bandwidth) interference is only efficient if the pathlengths 
are matched to within a wavelength or so - a stringent requirement. 



2.4-2. Field-of-View One consequence of the short coherency envelope for broadband 
observations is a limitation on the field-of-view. Bandwidth-smearing, as it is called, 
limits the field-of-view to be equal to the fringe-spacing x the number of fringes in 
the coherency envelope (see last subsection), FOV ~ Baseline ^ '^^dians. This "field- 
of-view" is thus baseline-dependent, leading to confusing interpretations of data for 
extended sources. While this effect can be modeled, it should be avoided by using a 
spectrometer to limit the bandwidth of individual observing channels. 

Another common limitation of the field-of-view is the primary beam of an individual 
telescope. For most kinds of beam-combiners, flux outside the diffraction-limited beam 
is rejected (spatial filtering is described more fully in §3.5.1). For most astronomical 
objects observed by interferometers, this is not a serious problem since long integrations 
by (low-resolution) individual telescopes can be used to confirm that no significant fiux 
arises from outside the primary beam. For an imaging interferometer, one would like to 
use narrow enough bandwidths so that bandwidth-smearing (on the longest baselines) 
is small enough so that the entire primary beam can be mapped. The requirement for 

this is approximately: ~ Longest Baseline ^ 

■"^■■^ AA ielescope Diameter 

Of course, having a wide field-of-view would be useful for many studies, such as 
measuring proper motions of stars at the galactic center. As discussed later, a wide field- 
of-view (beyond the primary beam) can only be achieved in a so-called Fizeau combiner. 
The Large Binocular Telescope Interferometer is the only interferometer currently being 
built which will have this unique and potentially very powerful faculty. 

2.4-3. Filling the (u,v) Plane The ability to make an image depends most strongly 
on the filled fraction of the (u,v) plane. Recall that the visibility amplitude and phase 
measured by an interferometer is directly related to a single component of the Fourier 
Transform of the object brightness distribution. If the object brightness is specified 
on coordinates of Right Ascension (pointing East) and Declination (pointing North), 
then the reciprocal Fourier space has axes referred to as (u,v). Following astronomical 
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notation, the positive u-axis typically points to the left on a diagram, just as right- 
ascension coordinates increase towards the left (East). 

For a fixed geometry of telescope locations, the Fourier coverage varies as the 
star rises and sets, and as a function of the star's declination and the interferometer's 
latitude. Figure 8 shows the Fourier coverage of three actual interferometers (number 
of telescopes 3, 6, and 21 for IOTA, CHARA, and Keck aperture masking respectively) 
for a declination 45° object spanning 3 hours before and after transit (assuming 
monochromatic light). In aperture masking, the pupil plane of a single telescope is split 
up into sub-pupils which are allowed to combine, just like a long-baseline interferometer. 
It is obvious that that the coverage increases rapidly with number of telescopes. It is 
not so obvious that the number of closure phases/triangles also rapidly increases with 
array size (equivalent to filling up the hyper- volume (ui,Vi,U2,V2) with closure triangles; 
see Table 1). Tuthill and Monnier (2000) studied how imaging fidelity and dynamic 
range are affected by differing amounts of Fourier coverage using real data. Obviously 
for imaging it is absolutely critical to collect as much coverage as possible, and suitable 
array design is further discussed in §3.2. 

3. Basic Designs of Stellar Interferometers 

3.1. Brief Historical Overview 

Here I give only a brief historical overview of progress in stellar interferometry drawn 
partially from the review by Lawson (2000a); I refer the interested reader to the above 
article for more information. 

Modern interferometry can be traced back to 19th century France. Hippolyte Fizeau 
first outlined in 1868 the basic concept of stellar interferometry how interference of light 
could be used to measure the sizes of stars. The first attempts to apply this technique, 
akin to modern-day "aperture masking," were carried out by E. Stephan soon thereafter, 
although the telescopes of the time had insufficient resolution to resolve even the largest 
stars. 

Albert Michclson developed a more complete mathematical framework for stellar 
interferometry in 1890; while apparently Michclson was unaware of Fizeau's earlier 
work, more historical investigation is needed to establish this definitively Along 
with Pease, Michclson (Michclson and Pease, 1921) eventually succeeded in measuring 
the diameter of a Orionis (Betelgeuse) in 1920-21 using the Mt. Wilson 100" 
telescope (following earlier measurements of Jupiter's moons; Michelson, 1890, 1891). 
Interestingly, Michelson needed a baseline longer than 100" in order to resolve Betelgeuse 
(uniform disk diameter ~47 mas) , and acquired one by installing a 20-foot interferometer 
beam on the Cassegrain cage as illustrated in Figure 9, reproduced here from their 
original paper. Following the success of the 20-foot interferometer. Pease (with Hale) 
constructed a 50-foot interferometer (on Mt. Wilson, but separate from the 100" 
telescope); although some results were reported, this experiment was not very successful. 



Optical Interferometry in Astronomy 



21 



An IGffA Configuration) 



CHARA Confijjrafon 






1 1 1 




O 1 


o 


o 




o 




o 



KeckAperUre Ma^ 




-am -im 



IOTA Fcurier Coverage 



CHARA Fourier Cow rg^e 




Keck Fourier Coverage 



' 



i" 1 • 



Figure 8. Example of (u,v) plane coverage for different interferometers. The top 
panels show the interferometer array configurations, while the bottom panels show 
the corresponding (u,v) plane coverage. For the IOTA and CHARA interferometers, 
I have assumed a source at 45° declination observed for three hours both before and 
after transit. The right-most panels show instantaneous "snapshot" coverage for an 
optimized 21-telescope array, a geometry actually used in the Keck aperture masking 
experiment (Tuthill ct al., 2000c). Note that the circles in the top plot are not to the 
same scale as the individual telescopes diameters but have been enlarged. 



Due to its generally outstanding atmospheric conditions, Mt. Wilson continued to be 
a choice site for interferometry projects, subsequently hosting the Mark III, ISI, and 
CHARA interferometers. 

Following the disappointing results from the 50-foot interferometer, it would be 
decades before significant developments inspired new activity in the optical arena. 
Meanwhile, advances in radar during World War 11 spurred rapid development of radio 
interferometry. We refer the reader to Thompson et al. (2001) for a discussion of the 
development of radio interferometry beginning with the first radio interferometer built 
by Ryle and Vonberg in 1946. 

The unexpected success of "intensity interferometry" would inspire a host of new 
projects. The basic principle behind the intensity interferometer was laid out in Hanbury 
Brown and Twiss (1956a), and describes how correlations of intensities (not electric 
fields) can be used to measure stellar diameters. First results were reported soon 
thereafter (Hanbury Brown and Twiss, 1956b), leading to the development of the 



Optical Interferometry in Astronomy 



22 




Figure 9. This diagram from Michelson and Pease (1921, Figure 1) illustrates how a 
20-foot interferometer beam was installed on the Mt. Wilson 100" telescope in order 
to create, for the first time, an interferometer capable of measuring the diameter of 
stars beyond the Sun. Figure reproduced by permission of the AAS. 

Narrabri Intensity Interferometer. With a 188 m longest baseline and blue-sensitivity, 
this project had a profound and lasting impact on the field of optical interferometry, 
measuring dozens of hot-star diameters (e.g.. Brown et al., 1967a,b; Davis et al., 1970; 
Hanbury Brown et al., 1970, 1974a). The small bandwidths attainable with Intensity 
Interferometry limited the technique to the brightest stars, and pushed the development 
of so-called "direct detection" schemes, where the light is combined before detection 
to allow large observing bandwidths. This group would go on to develop the SUSI 
interferometer. 

Dr. Charles Townes, inventor of the maser, began a novel interferometer project 
during this same time period at University of California at Berkeley. He used heterodyne 
receivers as in radio interferometry, but the local oscillators were CO2 lasers operating at 
frequencies of ~27 THz (or ~10 //m wavelength), orders of magnitude higher than radio 
or microwave oscillators. First experiments were performed using the twin McMath 
auxiliary telescopes (separation 5.5 m) at Kitt Peak, AZ; first fringes were obtained 
on the limb of Mercury in 1974 (Johnson et al., 1974) and on stars in 1976 (Sutton 
et al., 1977, 1978, 1979, 1982) where cool dust shells were detected around many late- 
type stars, heterodyne detection also suffers from bandwidth limitations (like Intensity 
Interferometry) as well as additional noise contribution from laser shot-noise, which 
becomes progressively worse at higher frequencies. The Townes group went on to develop 
the ISI interferometer on Mt. Wilson. 

Other mid-infrared efforts also are worthy of note. An independent project at 
Arizona using a kind of aperture masking on a single large telescope (using direct 
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detection) took place almost simultaneously with the Townes' experiments (McCarthy 
and Low, 1975; McCarthy et al., 1977, 1978). In France, Jean Gay and collaborators 
pursued long-baseline interferometry, both through heterodyne detection (e.g.. Gay and 
Journet, 1973; Assus et al., 1979) and later direct detection efforts (e.g., Rabbia et al., 
1990). 

Long-baseline interferometry on a star by directly combining the electric fields 
before photon detection ("direct detection") was first accomplished in 1974 by Labeyrie 
(1975), using a 12-m baseline. This continued the long history of interferometry 
innovation in Prance (starting from Fizeau), and many important experiments have 
followed. I note that "Speckle Interferometry" was first described by Labeyrie (1970b) 
and these ideas were also very influential to the field. However, I will largely limit this 
review to separate-element, or long-baseline, interferometry, and will omit comments on 
speckle. Following this 1974 demonstration in Nice, the project moved to the Plateau 
de Calern site and become known as the Interferometre a 2 Telescopes (I2T). The I2T 
made measurements in the visible (e.g., Blazit et al., 1977) and in the near- infrared (di 
Benedetto and Conti, 1983; di Benedetto, 1985). The Grand I2T (GI2T, Mourard et al., 
1994) began soon thereafter and was developed in parallel with the I2T on the same 
plateau, but with larger telescopes (1.5 m) and longer maximum baselines (up to 65 m). 

At a time when it was a struggle to simply get two telescope interferometers 
working, considering the array of telescopes needed for imaging was indeed far-fetched. 
Thus, imaging using optical interferometry began with aperture masking experiments 
on large single- aperture telescopes (in the tradition of Michclson). In aperture masking, 
a pattern of holes (size <ro in diameter) is cut in a plate and placed in the pupil plane 
of a large telescope. The interference pattern formed thus simulates one from an array 
of telescopes combined like a Young's multi-sht experiment. Baldwin et al. (1986b) and 
Haniff et al. (1987b) showed how aperture masking in the visible yield data identical to 
that expected for an imaging array, and produced images of binary stars using closure- 
phase imaging. This group, based at the University of Cambridge, England, would soon 
begin developing the COAST interferometer, which would succeed in producing the 
first image with an aperture synthesis optical array (Baldwin et al., 1996b). Infrared 
aperture masking at the Keck Telescope (Tuthill et al., 2000c) also grew out of work from 
this group in collaboration with the U.C. Berkeley ISI team, and excited unexpected 
imaging results from this work have led to much enthusiasm for developing infrared 
imaging capabilities for long-basehne interferometers such as CHARA and VLTI. 

There is one remaining important interferometer lineage to mention, one which 
led to the modern development of "fringe-tracking" interferometers such as the NPOI, 
PTI, and Keck Interferometers, as well as numerous experimental innovations. The 
Massachusetts Institute of Technology and the Naval Research Laboratory built and 
operated a series of prototype interferometers, named the Mark I, Mark II, and the 
Mark III. Shao and Staelin (1980) reported the first successful active fringe-tracking 
results, and this group has been most active at pushing the use of interferometers 
for precision astrometry. The Mark III was located on Mt. Wilson and was a fully 
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automated interferometer operating in the visible with basehnes up to 20 m (Shao et aL, 
1988). The high efficiency aUowed many astronomical programs to be carried out until 
it was shut down in about 1993, and is widely considered one of the most productive 
interferometers to date. Numerous articles were published covering areas of astrometry, 
angular diameters, precision binary orbits, and limb- darkening (e.g., Mozurkewich et al., 
1988; Hutter et al., 1989; Mozurkewich et al, 1991; Armstrong et al, 1992; Hummel 
et al., 1995; Quirrenbach et al., 1996). 

I should also mention the prototype interferometer IRMA (Infra-Red Michelson 
Array), built at University of Wyoming (Dyck et al., 1993). While this instrument did 
not operate for very long, those involved were largely responsible for initial success with 
infrared observing at IOTA and have had lasting impacts at a number of other currently 
operating U.S. facilities, including NPOI, PTI, and Keck Interferometers. 

This section was meant to introduce historical interferometers (ones no longer in 
operation) which have had a lasting impact on the field, and I have left descriptions of 
currently operating interferometers to §3.6. As discussed in the opening, this review is 
not meant to document all the important results from these first generation facilities, 
but rather to give appropriate historical background for understanding the current state- 
of-the-field. 

3.2. Overview of Interferometer Design 

Compare Young's two-sht experiment (Figure 1) to what you see in Figure 10. We see 
telescopes instead of slits and a beam combiner (with relay optics) instead of a screen for 
viewing the fringes. In a real interferometer, we must use delay lines to compensate for 
geometrical delay introduced by sidereal motion of a star across the sky; in this way, we 
"point" the interferometer at the target. In order to successfully interfere light together, 
each interferometer will have many subsystems, and in this review we will describe the 
state-of-the-art developments for the Telescopes, the Relay Optics, the Delay Lines, and 
the Beam Combination. 

Before discussing each of the critical subsystems, the importance of the physical 
placement of the telescopes for imaging will be discussed. Many of these issues are 
discussed in more detail by Mozurkewich (2000), and here we consider array design from 
the perspective of imaging, not for specialized purposes such as nulling or astrometry. 

If there were no practical constraints and telescopes could be placed optimally, one 
could consider many possibilities. Studies have been published considering distributions 
based on optimizing uniformity of (u,v) coverage using three-fold symmetric patterns 
(used in Keck aperture masking, Golay, 1971), Reuleaux triangles (used for the Sub- 
Millimeter Array, Keto, 1997), a spiral zoom array (considered for the Atacama Large 
MiUimeter Array, see ALMA memos #216, 260, 283, and 291), and a Y-shaped array 
(adopted by the Very Large Array). While the first three of these methods offer better 
Fourier coverage than the Y-shaped array, the "imaging" interferometers of NPOI and 
CHARA both use a Y-shaped array - why? 
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Figure 10. This schematic ihustrates the major subsystems of a modern optical 
interferometer: the Telescopes, the Relay Optics, the Delay Lines, and the Beam 
Combination. 



For the VLA, an important reason for using a a Y-shaped array was a practical one; 
it was easy to move the telescopes along railroad tracks in order to cheaply and easily 
reconfigure the array geometry. While optical telescopes in arrays do not generally run 
on tracks (except at IOTA), the desire to transport light to a central facility (see §3.3.2) 
leads one to a Y-shaped geometry where the three-arms of the array are defined by 
vacuum pipes which relay beams from the telescopes to the delay lines and combiners. 
The NPOI interferometer (shown here in Figure 11) has many stations along the three 
vacuum arms where telescopes can be located, thus creating a flexible, rcconfigurable 
system capable of pursuing many astronomical programmes. Another theoretical beneflt 
of this design is that many telescopes can be arranged along each arm allowing "baseline 
bootstrapping" for imaging highly resolved targets, a technique where strong fringes 
measured between close-by telescopes are used to "phase-up" the fringes on the longer 
basehnes. 

3.3. Critical Subsystems Technologies 

3.3.1. Telescopes All interferometers need light collectors of some kind. In many cases, 
simple "siderostats" are used, whereby a steerable fiat mirror directs starlight either 
directly to the interferometer or first through a beam-compressor ( "afocal" telescope) . 




Figure 11. Overhead view of the NPOI interferometer array, to show the Y-shaped 

array layout, defined by vacuum pipes extending out to many possible siderostat 
"pads," or stations. Photograph reproduced with permission of the Naval Research 
Laboratory. 



A siderostat has limited sky coverage and makes polarization measurement difficult 
(due to the changing, non-normal reflection angles off the flat), but is thought to offer 
a more stable structure for minimizing vibrations and pivot-point drifts for accurate 
astrometry. More recent interferometers, such as CHARA, VLTI, and Keck, have chosen 
traditional altitude- azimuth ("alt-az") telescope designs which give full-sky coverage and 
potentially salvaging polarization work. In addition, all interferometer telescopes have 
incorporated fast "tip-tilt" guiding which tracks (and corrects) fast jitter of the stellar 
image, usually using visible-light "quad-cell" detectors. This corrects the first-order 
term of the wavefront perturbations, aligning the wavefronts to allow for stable beam 
combination. 

Without high-order adaptive optics, there is little use for a telescope aperture much 
larger than the atmospheric coherence length tq (see §2.2). Hence, most telescopes 
in today's interferometers are small by "modern" (8-m class) telescope standards. 
Dedicated visible-light interferometers (e.g., NPOI, SUSI) have telescope apertures 
around 12-14 cm in diameter; near-infrared interferometers (e.g., PTI, IOTA, COAST) 
have apertures diameters around 45 cm. The recently-built CHARA interferometer 
includes 1 m apertures which can take advantage of excellent seeing conditions in the 
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infrared, and could benefit from adaptive optics correction; the Keck and VLT auxiliary 
telescopes were specified to be 1.8 m for similar reasons. However, interferometry is 
not just for "small" telescopes anymore, since the world's largest telescopes, the two 
Keck telescopes and also the four VLT telescopes, are now part of the new generation 
of optical interferometers. As of 2002, only the Keck Interferometer has observed using 
adaptive optics, although the VLT Interferometer will soon possess this capabihty. See 
Table 3 for a summary of telescope apertures of today's interferometers. 

3.3.2. Relay Optics, Delay Lines, & Metrology After being collected by the telescopes, 
the light must be directed to a central facility for beam combination. While it may 
seem trivial to set up a series of mirrors for this purpose, there are many subtle issues 
that must be addressed. Traub (1988) discussed how the geometry of the relay optics 
must not corrupt the relative polarization of the beams, due to differential phase shifts 
between the s- and p-wave reflections from the mirror surfaces for non-normal incidence. 
One must pay attention to the issues of mirror and window coatings as well as geometry. 

In addition, due to the long path lengths between the telescope and central 
beam combining facility, significant differential chromatic dispersion occurs if the 
light is propagating in air. In order to combine broad bandwidths, one must either 
transport the light through a vacuum or construct a dispersion compensator (Tango, 
1990; ten Brummclaar, 1995), whereby wedges of glass are inserted into the beam to 
compensate for air's index of refraction; a combination of partial vacuum plus dispersion 
compensation is also possible. The size of the mirrors in this optics chain is also 
important for limiting the effect of diffraction (Horton et al., 2001), and often also 
sets the field-of-view of the interferometer. Lastly, because of the many reflections, 
high reflectivity of the relay optics must be maintained to maximize throughput and 
sensitivity; a side-beneflt of evacuated relay optics is that the mirrors stay clean. 

Because of the Earth's rotation, the apparent position of an astronomical object 
is constantly changing. In order to track this sidereal motion, a movable delay line 
is needed to compensate for changing geometrical delay between wavefronts reaching 
any two telescopes. The diagram in Figure 10 shows this delay line as a right-angle 
retroreflector, although most interferometers do not actually use this geometry. The 
requirements on this system are amazingly stringent: nanometer-level precision moving 
at high speeds (> 1 cm/s) and over long distances (>100 m) - a dynamic range of 
> IQi"! 

By far the most popular architecture today for the moving delay line is based on 
the solution implemented by the Mark III interferometer (Shao et al., 1988; Colavita 
et al., 1991). The retroreflection is produced by focusing the incoming beam to a point 
coincident with a small flat mirror (attached to a piezo-electric stack), which reflects and 
re-collimated; a practical optical system is illustrated in Figure 12. This mirror system 
is mounted on a flexible stage which can be translated using a voice coil. Lastly, this 
whole stage is mounted on a wheeled-cart, which is driven on a rail by linear motors. 
This system has three nested feedback loops, driven by laser metrology: precise sub- 
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Figure 12. Diagram of tlic most standard delay line architecture used in optical 
interferometry, originally from the Mark III interferometer. Figure reproduced from 
Shao et al. (1988, Figure 4) with permission of ESO. 



wavelength control is maintained by the piezo-driven small mirror, when this mirror 
exceeds its normal operating range (~50//m) then offsets are given to the voice-coil 
stage, and so on. This basic architecture is in use at most interferometers built in the 
last 10 years; alternate delay lines include floating dihedrals mirrors on an air table 
(IOTA) and moving the beam combination table itself (GI2T). 

3.3.3. Beam Combination & Fringe Modulation Once the beams have been delivered 
to a central combination facility and have been properly delayed, there are many ways 
to actually do the interference. Here I discuss image-plane and pupil-plane combination 
(also know as "Fizeau" and "Michelson" combination, respectively), along with spatial 
versus temporal modulation of the fringes themselves. 

Figure 13 show these two different ways of detecting fringes in a two-element 
interferometer. In one case, an imaging system is used to flll the image plane with the 
equivalent of Young's fringes. As one moves along the image plane, there is a different 
relative delay between the interfering beams, and hence the modulation (fringes). In 
this scheme, there is no need to actively modulate the fringes; in fact, atmospheric 
turbulence will introduce relative delays and cause the fringe pattern to "slide" back and 
forth, smearing out the fringes on short time-scales if not stabilized. This combination 
scheme most closely follows the "two-slit" interferometer analogy developed in earlier 
sections. 

The second scheme, and currently the most common one, is pupil-plane 
combination, or "Michelson" -style combination. Interestingly, this method is named 
after Michelson, not because of his stellar interferometry work (which used "Fizeau" 
combination, see §3.1), but because of the interferometer used in the Michelson-Morley 
experiment. In this method, the wavefronts from the two coUimated telescope beams 
are overlapped on a 50/50 beamsplitter. Depending on the phase relationship of the 
waves, differing amounts of energy will be transmitted or reflected at the beamsplitter. 
Single-pixel detectors can then be used to measure the energy on both sides of the 
beamsplitter (the sum of which is conserved). The popular adoption of this method 
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Figure 13. Diagram of Image-plane and Pupil-plane beam combination techniques. 
The left panel shows image-plane, or Fizeau combination, where light from the two 
telescopes are brought together in an image plane to interfere, just like Young's two-slit 
experiment. The right panel shows how pupil-plane (or "Michelson" ) interferometry 
superimposes the two coUimated beams at a beam-splitter. By modulating the time 
delay of one beam with respect to the other (e.g., with the delay line), the interference 
can be modulated and fringes detected using single-pixel detectors. 

results largely from the signal-to-noise benefits of using single pixel detectors. In order 
to measure the amplitude of the coherence, a dither mirror (often in the delay line) 
sweeps through a linear pathlength difference of many wavelengths. The white-light 
fringe, or interferogram, can then be recorded, as long as the scanning takes place faster 
than an atmospheric coherence time. 

When dealing with an array of telescopes, there are more options, ten Brummelaar 
(1993) outlined some forward-looking beam combiner designs in the context of the 
CHARA array and useful articles by Mozurkewich (2000) and Mariotti et al. (1992) 
also contain extended discussion on the subject; here I only mention the highlights. 
The image-plane method can be extended to arbitrary number of telescopes, as long 
as the spacings between the beams are non-redundant, so that each beam-pair will 
have a unique fringe spatial frequency in the image-plane. Labeyrie (1996) elaborates 
on the concept of pupil densification, an idea finding increasing application in modern 
interferometry. The pupil-plane method can also be extended, either by combining 
the beams "pair- wise" or "all-in-one." In a pairwise-scheme, each telescope beam is 
split using beamsplitters and then various combinations are created to measure all the 
baselines. In the all-in-one scheme, more than two beams are superimposed and the 
fringes from different pairs are distinguished by modulating the delays such that each 
baseline pair has a unique fringe temporal frequency in the readout. There are methods 
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which combine pair-wise with all-in-one and are called "partial-pairwise." Each method 
has its advantages and disadvantages, depending on the availability of focal plane arrays, 
the level of readnoise vs. photon noise, the required calibration precision, etc. However, 
in general, "pair-wise" detection is the worst method for large number of telescope 
because the light has to be split more times (sec Buscher, 1988, although beware of 
some important simplifications made in this analysis). 

Coherent beam combination can be discussed more generally depending whether 
the interference occurs in the image/pupil plane and whether the telescope beams are 
co-axial or multi-axial. I refer the reader to the influential internal ESO report by 
Mariotti et al. (1992), which explains and defines the useful vocabulary in common use 
by the European interferometry community. 

We contrast the many imperfect beam combination strategies in the optical with 
those adopted in radio interferometry. At radio and microwave frequencies, the signals 
from each telescope can be split and re-amplified without introducing additional noise 
after the initial coherent detection (radio interferometers do not operate close to the 
Poisson limit). Hence, a pair- wise combination scheme can be employed without any 
loss in signal-to-noise ratio. In addition, the electric field at each telescope can be 
truly cross-correlated with that from all others leading to a kind of Fourier Transform 
spectroscopy. Further, this can all be done using digital electronics after fast digitization 
of the signals. For more information, see the description of the Hat Creek mm-wave 
correlator by Urry ct al. (1985). At the end of this process, the digital correlators 
can recover all baselines with arbitrary spectral resolution without lost sensitivity - a 
dramatically superior situation than possible in the photon-starved visible and infrared 
regime! 

3.3.4- Fringe Tracking An increasingly popular and powerful capability for optical 
interferometry is called "Fringe Tracking." To do so, the white-hght fringe has to be 
actively tracked because atmospheric fiuctuations cause the location of the fringe to 
vary by up to hundreds of microns on sub-second time scales. There are two levels 
of tracking these fringes, one is called "cohcrcncing" and the other is called "fringe 
tracking," although these terms arc often used rather loosely. 

In "coherencing," the interferometer control system will track the interferogram 
location to a precision of a few wavelengths. In a scanning interferometer, this will be 
sufficient to keep the full interferogram within the scanning range of the delay line. In 
an image-plane combiner, this will ensure you are near the peak of the white-light fringe 
(inside the coherence envelope set by the spectral bandpass). This can be done on a 
rather leisurely timescale, since large optical path distance (OPD) fluctuations tend to 
occur of slower time scale: update rates of ~1 Hz are sufficient except for the worst 
seeing conditions. 

True "Fringe Tracking" (also called "co-phasing") requires tracking OPD 
fiuctuations within a small fraction of wavelength in real-time, and hence requires 
orders of magnitude faster response (a time scale which depends on the wavelength and 
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seeing conditions). In the most common implementation (the "ABCD" method; see 
Shao and Staehn, 1977), two beams are combined pairwise while a mirror is stepped at 
quarter-wavelength intervals. The broadband white-light fringe is detected at one of the 
beamsplitter outputs, and fringe data is recorded synchronous with the dither mirror, 
resulting in four measurements (A,B,C,D) representing four different fringe phases. A 
discrete Fourier Transform (effectively) can be rapidly apphed to the data, resulting 
in a fringe phase estimate. This offset can be sent to the interferometer delay hne 
control system to nearly instantaneously correct for atmospheric turbulence (details 
in Shao et al., 1988; Colavita, 1999). The hght from the other beamsplitter output 
is usually dispersed and multi-wavelength data is collected. I also refer the reader to 
Lawson (2000), where the ABCD method (and other "phase" estimators) or compared 
to to "Group Delay" tracking methods, which use phase measurements at different 
wavelengths to measure interferometer delay offsets. 

Historically speaking, active fringe tracking has been important only for the 
Mark III interferometers and its successors (NPOI, PTI, Keck Interferometer). One 
reason fringe tracking has not been more widely pursued is because the sensitivity limit 
of a fringe tracking interferometer is less than fringe-envelope scanning interferometer. 
This is because very short integration times are required to stay on the fringe and hence 
the source must be fairly bright; in the fringe envelope scanning method, one has to only 
keep the interferogram in the scanning range and thus any given fringe measurement 
can have a lower signal-to-noise ratio. In practice, this amounts to sensitivity difference 
of a few magnitudes. 

As interferometers become more powerful and seek greater capabilities, fringe 
tracking is becoming a standard feature. High spectral resolution interferometry data 
is possible with fringe tracking systems, because a broadband white-light fringe can 
be used for fringe tracking while the remaining output channels can be dispersed. 
Normally, this data would have very low signal-to- noise ratio, but if the fringe tracking 
essentially "freezes" the turbulence, the dispersed fringes can be detected by integrating 
on the detector much longer than the typical atmospheric coherence time. Hence, fringe 
tracking is a kind of "adaptive optics" for interferometry. 

3.3.5. Detectors The most desired properties for detectors used in optical 
interferometry are low noise and high readout speed, two qualities usually not found 
at the same time. At the beginning of optical interferometry, the only visible-light 
detector was photographic film and infrared detectors were only just invented. Detectors 
have made incredible advances over the last few decades, and are operating near their 
fundamental limit in most wavelength regimes (the near- infrared is a notable exception). 

After years of struggling with custom-built photon-counting cameras for visible- 
light interferometry work, such as the PAPA camera (Papaliolios et al., 1985; Lawson, 
1994) and intensified CCDs (e.g., Blazit, 1987; Foy, 1988), commercial devices are being 
sold aimed at the adaptive optics market which have high quantum efficiencies (>50%), 
kilohertz frame times, and read noise of only a few electrons (fast readout CCDs). For 
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beam combinations schemes where single pixel detectors arc suitable, Avalanche Photo- 
Diodes (APDs) have as high quantum efficiency as CCDs but can photon-count at rates 
up to 10 MHz, although the best commercial devices seem to have an expensive tendency 
to stop working. This covers wavelengths from the blue to approximately the silicon 
cutoff (~1 iim). 

In the near-infrared (1-5/xm), there has been amazing progress this decade. 
After early work with single-element detectors (e.g, using material InSb), modern 
interferometers have taken advantage of technology development at Rockwell in near- 
IR focal plane arrays made of HgCdTe, such as the NICM0S3, PICNIC, and HAWAII 
chips. These arrays have high quantum efficiency (>70%) and can be clocked at ~MHz 
pixel rates with as low as 15 e- readnoisc. While not optimal, this represents orders-of- 
magnitudc improvement over photodiodes and has allowed new kinds of astronomical 
sources to be observed (most notably, young stellar objects). The noise can be further 
reduced by reading each pixel many times, a novel mode known as "non-destructive" 
readout. Hence, by reading a pixel n times before resetting, one can reduce the effective 
readnoise by approximately ^/n, for n<20. Interferometry benefits greatly from this 
capability, since only a few pixels need be readout, allowing large number of "reads" to 
be made in a short period of time (interferometry reference Millan-Gabet et al., 1999a). 

Traditionally, the HgCdTe detectors had a cutoff wavelength of 2.5^itn, but recent 
Molecular Beam Epitaxy (MBE) processes allow this cutoff to be tuned to much longer 
(or shorter) wavelengths (allegedly even beyond 5fim). For a 2.5fim cutoff, these 
detectors must be operated at liquid nitrogen temperatures (77 K) in order not to be 
saturated with dark current from thermally-generated electrons. An important recent 
development is that Raytheon has begun to compete with Rockwell in this market, and 
we can hope for even greater advances in HgCdTe arrays in the coming years as well as 
possibly even price reductions. 

Other materials, such as InSb, can be used for even longer wavelength performance. 
At 5/im and longer wavelengths, thermal background levels are sufficiently high that 
these detectors must be readout very rapidly, and can usually work in background- 
limited mode (despite >500 e- readnoise). This means that Poisson fluctuations in 
the thermal background flux dominate over other sources of noise (e.g., read noise); 
the only way to reduce the effect of this background noise is to increase the quantum 
efficiency of the detector or to reduce the thermal background load on the detector. 
Various companies have sold focal plane arrays in the "mid-infrared" (~8-25/xm) over 
the years, and are not all independent efforts after a complicated series of company sales 
(e.g., Hughes, Santa Barbara Research Center, Raytheon, Boeing). Recently Raj^heon 
has been offering Si:As Impurity Band Conduction (IBC) 320x280 focal plane arrays, 
which also operate at the background-limit. These detectors must be cooled below 
77 K to avoid high dark currents, and generally use liquid Helium. Uniquely, the ISI 
interferometer uses a single-element HgCdTe photodiode with high signal level (using 
CO2 laser local oscillator) which have up to 25% quantum efficiency and a 5 GHz output 
bandwidth. 
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3.3.6. System Control It is not trivial to control all the important subsystems of an 
interferometer. Many current interferometers (e.g., ISI, IOTA, PTI, Keck, VLTI) use 
the VME realtime architecture under the vxworks operating system (Wind Rivers). 
This allows different subsystems to be easily synchronized at the millisecond (or better) 
level. VME systems are fairly expensive, and some groups (in particular, CHARA) have 
adopted the RT (RealTime) Linux OS running on networked personal computers. 

3.4. Sensitivity 

Optical interferometers are orders-of-magnitude less sensitive than single-dish 
telescopes. At visible wavelengths where sensitivity is the worst, current interferometers 
have a similar limiting magnitude as the human eye (e.g., V mag ~6 at NPOI). In this 
section, we explore the current and future sensitivities of optical interferometers. 

3.4.1. What sets the limiting magnitude? There are three major problems which 
limit the sensitivity of today's interferometers: the atmosphere, optical transmission, 
detector/background noise. 

The sensitivity is most dramatically limited by the atmosphere which restrict the 
coherent aperture size and coherent integration time. We can use the notion of a coherent 
volume of photons which can be used for interferometry, with dimensions set entirely by 
the atmosphere. The coherent colume has dimensions of ro x ro x ctq, and hence is very 
sensitive to the seeing. Consider average seeing conditions in the visible (ro ~ 10 cm, 
to ~ 10 ms), we can estimate a limiting magnitude by requiring at least 10 photons 
to be in this coherent volume. Assuming a bandwidth of 100 nm, 10 photons (A ~ 
550 nm) in the above coherent volume corresponds to a V magnitude of 12.6, which 
is more than 10 magnitudes brighter than faint sources observed by today's 8-m class 
telescope. Because the atmospheric coherence length and time scale approximately like 
for Kolmogorov turbulence, the coherent volume oc A^. 

Current interferometers can not achieve this limiting magnitude because of 
additional problems. The most important in the visible is low optical throughput 
due to the large number of reflections between the telescopes and the final detector. 
The number of reflections easily exceeds 10 and is often closer to 20. Even with 
high quality coatings of 97% reflectivity, we see that ~50% of the light would be lost 
after 20 bounces (0.97^'^ = .54). In practice, current interferometers have visible-light 
transmission between 1% and 10%, due the fact that coatings degrade with time, the 
need for dichroics and filters with relatively high losses, and some diffractive losses 
during beam transport. Of course, detectors also do not have 100% quantum efficiency. 
The COAST interferometer has achieved the faintest limiting magnitude in the visible 
of ~9 mag, by optimizing throughput, detector quantum efficiency, and bandwidth 
(as another example, the NPOI interferometer which fringe-tracks and uses narrower 
bandwidths has a limiting magnitude around ~6). 

Throughput issues can be improved multiple ways. Lawrence Livermore Laboratory 
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is researching new coatings for mirrors which will have >99% reflectivity at most near- 
and mid-infrared wavelengths. In addition, simplified beam trains with few reflections 
are being designed for next generation interferometers. Lastly, the use of fiber and 
integrated optics could potentially lead to high throughput systems in the future; these 
developments are discussed more fully in §3.5. 

The last major hmitation is noise associated with the detection. Some visible 
light detectors, such as the photon-counting Avalanche Photo-Diodes, are almost 
perfect in this regard, boasting very low "dark counts" (<100 ct/s) and high quantum 
efficiency. However, this is not true in the infrared. Even the best infrared detectors 
have ~10 e~ noise per read. While normal (incoherent) astronomers can afford to 
integrate for minutes or hours to collect photons, interferometrists must readout pixels 
within the atmospheric coherence time and thus are strongly limited by readnoise. As 
one moves further into the infrared (5-10/im), then thermal background fiuctuations 
dominate the noise budget. Again, the relatively short coherence times of atmospheric 
turbulence directly result in a poor limiting magnitude compared to incoherent detection 
(i.e., photometry). The best published near- infrared performance of a two-element 
interferometer was reported by IOTA (MiUan-Gabet et al, 1999a) using a NICM0S3 
detector: J mag (1.25//m) 6.9, H mag (1.65//m) 6.9, and K' mag (2.2/im) 6.2, where 
J, H arc dominated by readnoise and K' is dominated by fiuctuations of the thermal 
background. Soon, these limiting magnitudes will be eclipsed by the adaptivc-optics- 
corrected Keck and VLT Interferometers which should be able to observe fainter than 
10th magnitude. 

There is not much experience yet with mid-infrared observations using direct 
detection. The ISI heterodyne interferometer has observed stars as faint as ~360 Jy 
(LkHQ;,101 Tuthill et al., 2002), corresponding to a N band mag of ~-2.2, limited largely 
by narrow bandwidths (AA ~ 0.002//m). The VLTI mid-IR instrument MIDI wiU be 
capable of broadband combination and is forecast to have a limiting magnitude of ~1 Jy 
(N band mag ~4) using the 8 m VLT telescopes (assuming the thermal background 
fiuctuations can be well-calibrated for systematic errors). Shortly before this article 
went to press, VLTI reported first fringes with the MIDI instrument. 

A number of new technologies are being explored to push down the limiting 
magnitude of optical interferometers, and some of these are described in the next section. 

3.5. New Technologies and Techniques 

One exciting aspect to the field of optical interferometry is the aggressive implementation 
of new technologies to extend the limits of the sensitivity and calibration precision. In 
this section, I will discuss new developments which are impacting optical interferometry. 

3.5.1. Spatial Filtering and Single-mode Fibers The idea to use single-mode fibers in 
optical interferometry was originated by Proehly (1982), and work began to implement 
these ideas in both Prance (e.g., Connes et al., 1987; Reynaud et al., 1992), and in the 
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Figure 14. This figure shows how the FLUOR beam combiner uses spatial filtering 
and photometric monitoring to allow precision calibration of fringe visibilities. Figure 
reproduced from Coude Du Foresto et al. (1997, Figurel) with permission of ESO. 



United States (Shaklan and Roddier, 1987; Shaklan, 1989). Following initial fringe 

detection in 1991 using the Kitt Peak McMath telescopes (Coude du Foresto and 
Ridgway, 1992), the FLUOR experiment as implemented on the IOTA interferometer 
was a real breakthrough; the amazing improvement in calibration precision was 
documented in Coude Du Foresto et al. (1997) and Perrin et al. (1998). Currently, the 
advantages of spatial filtering are being implemented at virtually all interferometers, 
and here I briefly explain why it is so important. 

Figure 14 shows a schematic of a fiber-based interferometer, as sketched by Coude 
Du Foresto et al. (1997). When coupling starhght into a single- mode fiber, the couphng 
efficiency depends on how coherent the wavefront is from an individual telescope 
(Shaklan and Roddier, 1988). A single mode fiber thus essentially converts phase errors 
across the telescope pupil into amplitude fluctuations in the flber. Once coupled into 
the single-mode flber, the light can be partially split in order to monitor the amount of 
coupled light ("photometric" outputs), and also can be interfered with light from another 
flber using a coupler, the flber equivalent of a beamsplitter. The transfer function of 
the fiber coupler is very stable and not dependent on the atmosphere; only the input 
coupling efficiency at each fiber is dependent on the atmosphere. Hence, the visibility 
can be measured very precisely (<0.4% uncertainty on reported by Perrin, 2003) by 
measuring the fringe amplitude and calibrating with the "photometric" signals. 

This method strongly mitigates the dominant source of calibration error in most 
optical interferometers, the changing atmosphere. The atmospheric turbulence must 
be monitored in some way when observing with an interferometer, since the coherence 
between two wavefronts from two telescopes strongly depends on seeing. However, this 
is not easy to measure with a typical interferometer in realtime, and hence one must 
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settle for interleaving "science" targets with "calibrator" sources to calibrate seeing 
drifts during the night. With fibers, the changing seeing conditions directly cause 
variations in the fiber coupling efficiencies which are monitored in realtime and corrected 
for. Figure 15 shows near-infrared visibility data on the calibrator star a Boo using 
both "conventional" interferometry and the FLUOR fiber optics beam combiner. The 
improvement to calibration is indeed dramatic and has had far-reaching effects on the 
direction of the whole field of optical interferometry. 

There are other ways to implement these calibration advantages than the FLUOR 
method shown in Figure 14. Monnier (2001) showed how the signal-to-noise can be 
somewhat improved by using an asymmetric coupler instead of separate photometric 
signals. Also, Keen et al. (2001) compared single-mode fibers with spatial filtering by 
small pinholes in order to determine which method is superior under different conditions. 
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Figure 15. (left panel) a. This figure shows visibility data for a Boo by the CERGA 
interferometer (diamonds) and the IRMA interferometer (squares), and originally 
appeared in the Publieations of the Astronomical Society of the Pacific (Copyright 
1993, Astronomical Society of the Pacific; Dyck et al., 1993, reproduced with permission 
of the Editors), b. The incredible gain in calibration using spatial filtering and 
photometric monitoring is evident in this figure reproduced from Perrin et al. (1998, 
Figure 2a) with permission from ESO. 



I should emphasize that there are many problems and limitations associated with 
using single-mode fibers, most notably low coupling efficiencies, high dispersion and poor 
polarization stability. Such problems have kept fiber optics from playing an important 
role for beam transport (Simohamed and Reynaud, 1997), and currently fibers are used 
only for beam combining and spatial filtering at specific wavelengths. For instance, 
sihca-based (telecom) fibers can generally only be used at J (1.25/im) and H (1.65/im) 
bands; the FLUOR experiment utilized Fluoride glass fibers which can transmit at K 
band (2.2^m) and beyond. Advances in the field of photonic crystals and photonic 
bandgap materials could lead to new fibers with low dispersion and high transmission 
for new interferometry applications, and should be aggressively pursued. 
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3.5.2. Integrated Optics While combining two telescopes together is straightforward 
using fiber optics, it becomes very difficult for multiple telescopes. This is because the 
light has be split many times, combined together many times, and the fiber lengths must 
be precisely matched and maintained to correct for differential chromatic dispersion and 
birefringence effects. 

An elegant solution to this problem, while maintaining the advantages of spatial 
filtering, is the use of integrated optics, the photonics analog to integrated circuits. 
P. Kern and an active group centered at Grenoble Observatory have pioneered this 
technique (e.g.. Kern et al., 1997; Malbet et al., 1999; Berger et al., 1999) and it is finding 
successful application at a number of observatories, including IOTA (Berger et al., 2001) 
and VLTI. In these combiners, many fibers can be mated to a small planar element 
with miniature waveguides etched in place to manipulate the light (split, combine, etc). 
Dozens of beamsplitting and combinations can all be fit into a few square centimeters 
- and never needs re-aligned! 

While integrated optics can solve the problem of how to combine many beams 
using guided optics, it has similar difficulties as fibers of poor transmission, limited 
wavelength coverage, dispersion, and birefringence. While the commercial applications 
for integrated optics in telecommunications has driven much of the innovation in this 
field, the astronomy community must actively engage with the photonics engineers 
to design custom components which can overcome the remaining problems for next- 
generation "astronomical-grade" devices. 

3.5.3. Adaptive Optics One critical advance to improve the sensitivity of infrared 
interferometers is the application of adaptive optics on large aperture telescopes. 
Generally, visible light photons are used to measure the wavefront distortions in realtime, 
allowing them to be corrected using a deformable mirror. Once the aperture is "phased- 
up," the entire (much larger!) coherent volume can be used for the infrared fringe 
detection. This method has already been applied on the Keck Interferometer, where 
AO systems on the individual 10 m telescopes now allow observations approaching K 
mag 10 (and should allow even fainter objects eventually). The major drawback for this 
is that there must be a "bright" visible guide star in the isoplanatic patch for the AO 
system to use for wavefront sensing, not possible for obscured sources such as Young 
Stellar Objects and dusty evolved stars where the visible source is often too faint (a few 
AO systems do have infrared wavefront sensors to mitigate this problem, e.g. Brandner 
et al., 2002). The maturation of laser guide star adaptive optics will allow this gain 
in coherent volume for all infrared observing eventually. Of course, building future 
interferometers at the most excellent sites (even in space) will be an increasing priority. 

3.5.4- Phase Referencing Phase Referencing is a kind of adaptive optics for 
interferometry, where a bright reference star is used to measure and correct for 
atmospheric time delays. This technique is used in radio interferometry to allow long 
coherent integrations on targets, via a fast switching scheme. 
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Unfortunately, the short atmospheric coherence times make a switching scheme 
difficult to implement. A different approach pursued by the Palomar Tcstbcd 
Interferometer is to use a "dual star module" , where light from two stars are selected and 
observed simultaneously using different delay lines. This allows both relative astrometry 
and phase referencing to be achieved. 

Very few results have been published on this technique so far, although the 
technique is being implemented at the Keck Interferometer and is planned for VLTI. 
First results from PTI have been pubhshed (Lane and Colavita, 1999; Lane et al., 2000a; 
Lane and Colavita, 2003), reporting extending the atmospheric coherence time to 250 ms 
and visibility calibration precision of 3-7%. Development of this technique will allow 
very faint limiting magnitudes, for a small set of sources located within an isoplanatic 
patch (~30" ) of a bright star. 

Another method called "differential phase" is being applied soon, where fringes 
at one wavelength are basically used to stabilize fringes at all the others. When a 
source shows significant wavelength-dependent structure, this technique should prove 
very powerful. This is discussed further in the context of extrasolar planet detection in 
§5.L3. 

3. 5. 5. Spectroscopy Very little has been done in the area of interferometric observations 
on spectral lines. The best science results will be reviewed in the next section, 
however here 1 wanted to mention recent developments. Bedding et al. (1994) discussed 
methods for combining aperture masking with spectroscopy, and the design of the 
MAPPIT instrument offers lessons for long-baseline interferometrists. G. Weigelt and 
collaborators have developed a spectrometer for use on two element interferometers, 
which allows near-infrared molecular bandheads of CO and H2O to be spatially probed 
(e.g. Weigelt et al., 2000; Hofmann, 2002). More interestingly, the AMBER instrument 
for the VLTI will boast three different spectral resolutions (up to R ~ 10000 across the 
infrared), making observations of individual lines possible (e.g., Petrov et al., 2000). 
While the G12T has had high spectral resolution for years, a number of other visible- 
light interferometers, including NPOI and COAST, have modified their combiners to 
allow Ha interferometry, following the fascinating results of the GI2T in the 1990s (e.g., 
Vakili et al., 1998); see §4.2.1 for more discussion on this. 

3.5.6. New Detectors Single-pixel visible light detectors are nearly ideal in their 
performance (e.g., APDs). However, new detectors exist with many of the same 
advantages of APDs, but which can also measure the energy of each detected photon 
(Superconducting Tunnel Junction detectors. Peacock et al., 1997). Although limited to 
maximum count rates of ~10 KHz, current STJ devices offer high quantum efficiency, 
timing accuracy, and about ~12% bandwidth energy resolution in the visible and have 
been used on the sky (Perryman et al., 2001). One obstacle for this technology is 
that most astronomers want large-format focal plane arrays with millions of pixels, and 
present arrays are ~6x6 pixels. These STJ arrays are small, but large enough to be be 
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quite interesting for optical interferometry; this work should be strongly encouraged. 

For some type of interferometer combiners (e.g., high-resolution spectrometers or 
6-telescope imaging), many pixels are needed; unfortunately, APDs and STJs are not 
economical for this and CCDs typically still have larger readnoise for fast frame rates. 
In this regard, a new development by Marconi may be interesting (Mackay et al., 2001). 
They have produced a kind of "photon-counting" CCD, which implements on-chip 
avalanche gain stages in order to amphfy single electrons into large signals. Tubbs 
et al. (2002) report the first use of these new detectors in astronomy, and the results 
are promising for interferometry (the COAST interferometer is currently adapting such 
a device for their work) . 

Because of the relatively high readnoise for near-infrared detectors, improvements 
in the next decade could easily extend the sensitivity of interferometers by a factor of 
10. The AOMUX detector program by Rockwell has just begun, and has the goal of a 
few electron readnoise at high frame rates. Keeping pace with these developments will 
remain a high priority for optical interferometry. 

There are also some developments to create photon-counting near-infrared 
detectors, equivalent to Avalanche Photo-Diodes. Sometimes called SSPMs (Solid 
State Photo-Multiphers), Eikenberry et al. (1996) described one experiment. Currently, 
the main drawback with these devices is the low quantum efficiency, a few percent. 
Alternatively, Superconducting Tunnel Junctions can also be used in the near-infrared 
for photon-counting. 

3.5.7. Nulling Another interferometric technique gaining application is nulling. By 
introducing an achromatic 180° phase shift in one beam, the white-light fringe can 
be turned into a white-hght null (Bracewell, 1978). This has obvious apphcations for 
extra-solar planet searches and zodiacal dust disk characterizations, since removing the 
bright central star is essential for detecting faint circumstellar material and companions. 
The only astronomical results from nulling have come from aperture masking style 
experiments (e.g., Hinz et al., 1998; Hinz, 2001; Hinz et al., 2001b), and have encouraged 
aggressive follow-up experiments. In particular, the Keck Interferometer is pursuing 
a mid-infrared nulling project (Scrabyn and Colavita, 2001) and nulling is a central 
operational mode for the Large Binocular Telescope Interferometer (Hinz et al., 2001a). 

3.6. Current and Future Facilities 

In Tables 2 & 3, I have summarized all the current and planned facilities (ground- 
based). Further discussion of the current field, including documentation of the rising 
trend of pubhcations, can be found in Ridgway (2000) where I have found some of the 
information for these summary tables. We note that links to all these interferometers can 
be found on the well-established "Optical Long-Baseline Interferometry News" website, 
maintained by Peter Lawson at NASA-JPL (http://olbin.jpl.nasa.gov). 

Each of the currently operating interferometers have unique capabilities and 
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achievements of note. The GI2T and ISI Interferometers are the longest operating 
interferometers, both beginning work in the 1980s; notably, the G12T has uniquely 
pursued observing of Ha emission (and remains the only direct detection interferometer 
with general high spectral resolution capabilities) and the ISI is the only (published) 
mid-infrared interferometer. The COAST and NPOI interferometers are currently best 
optimized for imaging, having incorporated 5 and 6 telescopes respectively into their 
arrays. IOTA is noted for groundbreaking fiber optics and detector development in 
the infrared. NPOI and PTI have incorporated elaborate internal metrology to enable 
ambitious astrometry goals. SUSI has the capability of 640 m baselines and is one of 
the only interferometers in the southern hemisphere. Strong progress from the MIRA-I 
array marks Japan's recent efforts in long baseline interferometry. 

Recent developments include new infrared and visible combiners for the IOTA (first 
integrated optics success with stars Berger et al., 2001) and GI2T interferometers, third 
telescope upgrade projects for the ISI and IOTA interferometers, 6-telescope operation 
by NPOI and 5-telescopes for COAST, and first fringes from the Keck, VLTI, MIRA, 
and CHARA interferometers. SUSI has also commissioned a new "red" table, allowing 
packet-scanning interferometry using APDs. The FLUOR combiner, so successfully used 
on the IOTA interferometer, has been moved to CHARA, and we can expect excellent 
results soon to take advantage of the greater resolution and sensitivity. 

Indeed, it has been a busy decade for construction and implementation. It 
is apparent in Table 3 that the current and next generation interferometers boast 
significantly larger and more numerous telescope apertures and baselines, and promise 
to deliver significant new results. In the area of imaging, CHARA and NPOI will have 6 
telescopes spread over hundreds of meters, to allow imaging capabilities at milliarcsecond 
resolution. The VLTI and Keck Interferometers will have >100 m basehnes with adaptive 
optics corrected primary mirrors, allowing many new kinds of science to be pursued. In 
particular, we can expect the first extragalactic sources, bright AGN and quasars, to be 
measured at near-infrared wavelengths very soon (probably before this article goes to 
press). These new developments are further discussed in §5. 

Lastly, I will mention recent progress on the next generation of interferometers. 
The GHANA project has carried out initial experiments to couple light from Mauna Kca 
telescopes into single-mode fibers, the first step in a plan to link the giant telescopes of 
the Hawaii into a powerful optical interferometer. Major construction for the Large 
Binocular Telescope (and Interferometer) has been progressing for many years and 
is in an advanced stage now. Importantly, the final design plans for the Magdalena 
Ridge Observatory are shaping up and site work for the ~10 telescope optical array is 
expected to begin soon. You can find more information on these ambitious projects in 
the interferometer summary tables as well. 

The next section will review the currently exciting results from optical 
interferometry, and give some indication of how the new facilities will impact many 
areas of astrophysics. 



Table 2. Current and Future Optical Interferometers: Basics (* indicates "in planning") 



Acronym 


riiii Aanic 


Lead lust it ut iou(s) 


Loeat iou 


St art 


CHARA 


Center lor Hign Angular Resolution Astronomy 


Georgia State University 


T\ T 1 T T T * 1 /~i A T T O A 

Mt. Wilson, CA, USA 


2000 


COAST 


Cambridge Optical Aperture Synthesis ielescope 


Cambridge University 


Cambridge, England 


1992 


GET 


Grand Interferometre a 2 Telescopes 


Observatoire Cote D'Azur 


Plateau de Calern, France 


1985 


IOTA 


Intrared-Optical Telescope Array 


Smithsonian Astrophysical Observatory, 


H Ti TT 1 • A n TTn A 

Mt. Hopkins, AZ, USA 


1993 






Univ. oi JVlassachusetts (Amherst) 






TOT 


T £ JO J.*lTj.r 1. 

mrrarea bpatial interrerometer 


TT • £ IT • J. Tl 1 1 

Univ. of Cahiorma at Berkeley 


"\ /Tj. THI'l A TTO A 

Mt. Wilson, CA, USA 


1988 


Keck-I 


Keck interrerometer (Keck-1 to Keck-11) 


NASA-JPL 


H T TT" TTT T TO A 

Mauna Kea, HI, USA 


2001 


1\ /TTT) A T 

MiKA-1 


Mitakc inirared Array 


National Astronomical Observatory, Japan 


Mitaka Campus, iokyo, Japan 


1 c\c\o 

1998 




Wavy Prototype Optical Interierometer 


Naval Research Laboratory, 


Til ^„j.„xr A T TTC A 

li'lagstan. A/, USA 


1994 






TT O "VT 1 /^l J_ 

U.S. Naval Observatory 






PTI 


Palomar Test bed Interferometer 


"NT A O A TT~»T 

NASA-JPL 


A yTj- T~» 1 /~1 A TTO A 

Mt. Palomar, CA, USA 


1996 


CITTCIT 
oUol 


Sydney University Stellar Interferometer 


Sydney University 


Narrabri, Australia 




VLTI-UT 


VLT Interferometer (Unit Telescopes) 


European Southern Observatory 


Paranal, Chile 


2001 


Keck* 


Keck Auxiliary Telescope Array 


NASA-JPL 


Mauna Kea, HI, USA 


~2004? 


LBTI* 


Large Binocular Telescope Interferometer 


LBT Consortium 


Mt. Graham, AZ, USA 


-2006 


MRO* 


Magdalena Ridge Observatory 


Consortium of New Mexico Institutions, 


Magdalena Ridge, NM, USA 


~2007 






Cambridge University 






OHANA* 


Optical Hawaiian Array for Nanoradian Astronomy 


Consortium (mostly French Institutions, 


Mauna Kea, HI, USA 


~2006 






Manna Kea Observatories, others) 






VLTI-AT* 


VLT Interferometer (Auxiliary Telescopes) 


European Southern Observatory 


Paranal, Chile 


~2004 
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Table 3. Current and Future Optical Interferometers: Capabilities (* indicates 
planned capabilities) 





Telescope 


Maximum 


Wavelength 


Acronym 


Number 


Size (m) 


Baseline (m) 


Coverage 


CHARA 


6* 


1.0 


330 


visible*, near-IR 


COAST 


5 


0.40 


47 (100*) 


visible & near-IR 


GI2T 


2 


1.52 


65 


visible, near-IR 


IOTA 


3 


0.45 


38 


visible, near-IR, Afim 


ISI 


3* 


1.65 


85 (> 100*) 


mid-IR 


Keck-I 


2 


10.0 


85 


near-IR, mid-IR* 


MIRA-I 


2 


0.25 


30 


visible 


NPOI 


6 


0.12 


64 (> 250*) 


visible 


PTI 


3 


0.40 


110 


near-IR 


SUSI 


2 


0.14 


64 (640*) 


visible 


VLTI-UT 


4 


8.0 


130 


near-IR, mid-IR 


Keck* 


4* 


1.8 


140* ? 


near-IR 


LBTI* 


2* 


8.4 


23* 


near-IR, mid-IR 


MRO* 


-10 


-1.5 


-1000 


visible, near-IR 


OHANA* 


~6 


3.5-10 


-1000 


near-IR 


VLTI-AT* 


3* 


1.8 


202 


near-IR, mid-IR 



4. Summary of Major Scientific Results 

This section is divided up into two major areas: astrophysics of stars and of circumstellar 
environments. Optical interferometers have made substantial contributions in each, and 
I will outline recent progress. 

4.1. Stellar Astrophysics 

Optical interferometry has made the greatest impact in the area stellar astrophysics, 
in particular the study of nearby single stars. This is not surprising, given the limited 
nature of single-baseline interferometers and the limited sensitivity of first-generation 
instruments. In the last decades of work, an impressive diversity of investigations have 
been carried out and here we document the most successful work. 

4.1.1. Stellar Diameters and Effective Temperatures One of the earliest identified 
applications for optical interferometry was directly measuring the effective temperature 
scale of stars. The effective temperature is defined such that 

L = AnaR'T^^ (8) 

where a is the Stephan-Boltzman constant, R is the radius of the star, and L is the total 
bolometric luminosity. Hence, by measuring the angular size of a star and the apparent 
luminosity, the effective temperature can be directly calculated. The above equation 
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is often rearranged in terms of directly observable quantities (independent of distance 
estimate) : 

refr = 2341(Fboi/^^)'/' (9) 

where Fboi is the total bolometric flux (in 10~^ ergs cm~'^ s~^) and Or is the angular 
diameter in milliarcseconds. Empirical calibration of the effective temperature as a 
function of spectral type is important since Teff is considered a fundamental parameter 
of a star, appearing throughout stellar astrophysics most notably on the Hertzsprung- 
Russell diagram. 

The survey of stellar diameters using intensity interferometry by Hanbury Brown 
et al. (1974a) still serves as the best resource for the effective temperature scale of hot 
main sequence stars. The technique of lunar occultations has traditionally been the 
other main method for high-resolution measurements of stellar sizes, as represented by 
the classic paper by Ridgway et al. (1980). 

There are now more than a hundred interferometer diameter measurements, and 
this progress is marked in Figure 16. Here we see one of the flrst major results 
from Michelson interferometry from the I2T/CERGA interferometer (di Benedetto 
and Rabbia, 1987), the effective temperature scale of giants. Next to it, is a more 
recent version of the same diagram showing the increase in the number of diameter 
measurements, compiled by van Belle et al. (1999). The effective temperature scale for 
late-type stars is now well-established, and available diameter data has been generated 
by many interferometers (e.g., Mozurkewich et al., 1991; Dyck et al., 1996; Perrin et al., 
1998; Nordgren et al., 1999). A recent cross-comparison found the datasets from different 
groups to be statistically consistent (Nordgren et al., 2001); but for the latest spectral 
types, the visible photosphere is affected by TiO absorption and it is believed that 
infrared sizes are more representative of the "true" photospheric extent (e.g., Dyck and 
Nordgren, 2002). 

While giant stars have made easy targets for interferometers due to their large 
angular sizes and high luminosities, the census of lower-mass dwarf stars and hotter 
main sequence stars remain incomplete. The PTI interferometer recently made first 
contributions to the study of K- and M-dwarfs by measuring the diameters of 5 
such stars (Lane et al., 2001), although greater precision (<1% diameter errors) is 
needed to stringently test theoretical models. The first scientific result from the 
VLTI interferometer recently contributed to this precious, limited dataset of M-dwarf 
diameters (Segransan et al., 2003), also finding sizes consistent with theory but lacking 
precision. To date, there are no published diameters with baselines as long as the 
Narrabri Intensity Interferometer (188m), although the CHARA interferometer did 
record fringes on a 330m baseline in 2001. The new long baseline capabilities of CHARA, 
NPOI, and SUSI should allow progress in some of these areas in the near future. 

4.1.2. Limb-Darkening, Atmospheric Structure Measuring a so-called "uniform-disk 
(UD)" diameter only requires a single visibility data point for an isolated star, by fitting 
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Figure 16. (left panel) a. This figure shows one of the first major Michelson 

interferometer results (by the CERGA/I2T interferometer), the effective temperature 
relations for giants by di Benedetto and Rabbia (1987, see their Figure 1), reproduced 
here with permission of ESO. (right panel) b. Here is reproduced Figure 2a from van 
Belle et al. (1999), with permission from the AAS, showing the huge increase in the 
number of diameter measurements that now can be used for empirically determining 
the effective temperature scale of giants. The X-axis label "V-K" refers to the 
brightness of the star at V-band (Ao=0.55/im) compared to K-band (Ao=2.2/im); 
redder stars have larger V-K color. 

a one-parameter model. However, stellar photospheres are known to limb-darkened 
due to optical depth effects, and thus UD diameters must be corrected to yield the 
correct physical photospheric size. The first attempts to directly measure this was done 
using the intensity interferometer on the AlV star Sirius A (Hanbury Brown et al., 
1974b), however the results suffered from large errors and were not very definitive. 
Other attempts have been made based on looking for wavelength-dependent angular 
diameters, a sign of limb-darkening (e.g., Ridgway et al., 1982; Mozurkewich et al., 
1991). 

The visibility curve of a uniform-disk star is related to the first Besscl function, 
and contains an ever decreasing scries of lobes, separated by nulls, as one observes with 
increasing angular resolution (see Figure 15 for a plot of the first two lobes of this curve). 
The main difficulty for limb-darkening studies is that the first "lobe" of the visibihty 
pattern for a star is insensitive to limb-darkening effects (mathematically, it probes only 
the 2nd moment of the brightness distribution; see Lachaume, 2003); measurements 
beyond for the first null must be made to unambiguously detect limb-darkening effects 
(see Figure 17 for representative visibility curves). However, the fringe contrasts at high 
spatial resolution are necessarily low, and thus it has been difficult to measure these 
effects. Precise measurements of Arcturus {a Boo, KlIII) were made by Quirrenbach 
et al. (1996) using the Mark HI employing a novel phase-referencing method to increase 
the signal-to-noise near the visibility null; they found reasonable agreement with model 
expectations. 

The NPOI interferometer succeeded the Mark HI, and expanded the wavelength 
phase-referencing techniques (using strong fringes at one wavelength to allow coherent 
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integrations on weak fringes). Hajian et al. (1998) used these advantages first for limb- 
darkening studies, followed more recently by Wittkowski et al. (2001). The latter 
paper presented data with spatial frequencies sampled well before and after the null 
and these exemplary results are reprinted here in Figure 17. Again, the atmospheric 
models were found to be in reasonable agreement with the interferometry results. A new 
generation of precision tests of stellar atmospheres are now being pursued using detailed 
radiative transfer modeling coupled with thoughtful interferometer measurements at 
specific wavelengths and with specific baseline coverage (e.g., Aufdenberg et al., 2002). 



Figure 17. This figure shows recent results from the NPOI interferometer, which has 
been optimized to investigate effects of limb-darkening on stellar profiles. A portion of 
Figure 3 from Wittkowski et al. (2001) has been reproduced here, with permission from 
EDP Sciences, showing that models accurately predict the amount of limb- darkening 
observed in this K5 star. 



Until recently, most interferometers had only 2-elements and had difficulty to 
directly measure the two-dimensional shape of stellar photospheres. In most cases, 
one assumes the star is circularly symmetric in order to interpret visibility data taken 
at differently projected baselines. The Palomar Testbed Interferometer made the 
serendipitous, although not unexpected in retrospect (Hanbury Brown et al., 1967), 
discovery that that the rapidly rotating star Altair is not circular. Figure 18 shows the 
oblate spheroid model of this source developed by van Belle et al. (2001). The oblatcness 
is caused by centrifugal "bulging" along the equator and these measurements offer an 
independent measure of the projected stellar rotational velocity v sin i. While red giants 
and supergiants had been known to deviate from circular symmetric (sec §4.1.5), this 
is the first main sequence star found to be non-circular; future "imaging" work should 
allow new probes of other rotational effects, such as gravity-darkening. 

In addition to standard limb-darkening profiles, one has to be concerned with the 
effects of molecular lines formed in the photospheres of cool giants (mostly M spectral 
types), especially for the Mira variables. It has been known for quite some time (first 
from speckle interferometry) that evolved stars appear much larger when observed in 
narrow spectral channels coincident with deep TiO bands in the visible regime (e.g.. 
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Figure 18. Oblate Spheroid Model for the photosphere of Altair, based on PTI data 
(from Figure 6 of van Belle et al., 2001), reproduced here with permission of the AAS. 



Labeyrie et al., 1977). The origin of this extension is obvious: at the observing 
wavelength, the optical depth unity surface is at greater distance from the star, and 
hence the apparent size is noticeably larger. Quirrenbach et al. (1993b) made the first 
systematic study of this effect for late-type stars covering a range of spectral types, and 
these results are shown here in Figure 19a. 

A more recent, and unexpected, discovery in this vein is that molecules with 
transitions in the near-infrared are causing large increases in the apparent sizes also. 
The effect of (most likely) unappreciated H2O lines in the coolest M-stars was uncovered 
by many groups using different approaches at about the same time. Perrin et al. (1999) 
detected hints of these effects, finding puzzling deviations from a uniform disk for the 
0-rich Mira R Leo. Tuthill et al. (2000a) found that R Aqr, another 0-rich Mira, was 
dramatically larger at 3.1 /im than at shorter wavelengths (an effect seen also in other 
0-rich Miras, Tuthill et al., 1999c). First results with an L' band (3.75/xm) combiner at 
IOTA also found a large diameter increase compared to shorter wavelengths (Mennesson 
et al., 1999). 

A possible explanation for this effect was separately noted by researchers analyzing 
data from the Infrared Space Observatory (ISO), finding new water features in this part 
of the spectrum (e.g., Tsuji et al., 1997; Matsuura et al., 2002). Another recent analysis 
(Jacob and Scholz, 2002) coupled a dynamical model to a simple radiative transfer 
model and found complex (time-variable) visibility curves due to molecular effects in 
the near-IR. Mennesson et al. (2002) have collected data from multiple infrared bands 
(e.g., Chagnon et al., 2002), arguing the presence of "extended gaseous layers" around 
0-rich miras; see the dramatic difference in near-infrared sizes observed for R Aqr in 
the right panel of Figure 19. 

Thompson et al. (2002) have expanded these studies by measuring the sizes of 0-rich 
and C-rich miras using narrow spectral channels (AA ~ O.l/xm) from 2.0-2.4/im. While 
they report 0-rich miras are larger near the edges of the band, a different behavior is 
observed for C-rich stars; this tentatively confirms the role of 0-bearing molecules (e.g.. 
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H2O). Moving further into the infrared, Weiner et aL (2000) (following earlier work by 
Bester et al., 1996) actually find stars to be larger in true 11.15/xm continuum channels 
(using the ISI interferometer) than at 2.2fim, a somewhat confusing result since the near- 
IR wavelengths are expected to be significantly contaminated by molecular effects. I am 
aware of even more results which have not made it to press yet (e.g., adaptive optics at 
the Subaru Telescope, narrow-band IOTA interferometry) , and anticipate rapid progress 
in this area over the coming years. 



Figure 19. (left panel) a. This figure shows the Mark III compendium of results 
measuring photospheric diameters in and out of strong TiO bands (reproduced from 
Figure 2 of Quirrenbach et al., 1993b, with permission of the AAS). Redder stars 
(later spectral types) show greater atmospheric extensions in TiO. (right panel) h. 
This figure shows the most recent data illustrating that Mira stars have strongly 
wavelength-dependent diameters (reproduced from Figure 2 of Mennesson et al., 2002, 
with permission of the AAS). The two curves were fit to visibility data taken at 2.2/xm 
and at ^4/im, and show a greater than factor of two increase in angular size between 
these two wavelengths. 



4.1.3. Pulsating Stellar Atmospheres As discussed earlier, the "continuum" 
photospheric size is important for calculating a meaningful effective temperature scale. 
Another important consequence of angular diameter work is specific to variable stars: 
the average physical size should reveal whether a star is pulsating in the fundamental or 
first-overtone mode. Distance estimates have been combined with interferometry data 
to estimate physical diameters; these studies typically found "large" sizes consistent 
with first-overtone pulsation in most cases (e.g., Haniff et al., 1995; van Leeuwen et al., 
1997; Whitelock and Feast, 2000), although some sources were found to be fundamental 
pulsators. This is at odds with both non-linear pulsation models of Miras (Bessell ct al., 
1996) as well as the persuasive study of variables in the Large Magellanic Cloud by 
Wood and Sebo (1996). If visible and near-infrared diameters are indeed contaminated 
by molecular absorption as indicated by recent intcrfcromctric results discussed in the 
last section, it is possible that the true continuum diameters are small enough to be 
consistent with fundamental mode pulsation. The pulsation mode question of Miras has 
been debated and "settled" many times, and still more work is needed for a definitive 
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answer. 

Pulsating stars, especially the Mira variables (period ~ 1 year), are also expected to 
have large changes in the photospheric diameters as the luminosity varies (e.g., Bessell 
et al., 1989; Ya'Ari and Tuchman, 1996). However, because of the difficulty in obtaining 
a uniform dataset over many years, it is only recently that good pulsation curves of 
diameters have become available. 

While first hints of phase-dependent diameter changes were reported based on a 
statistical analysis of IOTA data (van Belle et al., 1996), the first definitive detection 
of diameter pulsation came from the COAST group for the 0-rich mira R Leo (Burns 
et al., 1998), where a 35% change in diameter was reported. In Figure 20a, we show 
more recent results for x Cyg (also from COAST) with better temporal sampling of the 
pulsation curve (Young et al., 2000b). Most recently, Weiner et al. (2003) report the 
first detection of pulsation at mid-infrared wavelengths, in this case for Mira variable 
Ceti. 

The most significant recent developments are coming from the Ph.D. dissertation 
of Thompson at PTI. The high level of automation of the PTI has allowed systematic 
observations of a large number of Miras at all pulsational phases, and these data are 
presented in Thompson (2002). A first look at the data has been pubhshed, and one 
result from Thompson et al. (2002) appears in Figure 20b; the diameter of this star has 
been measured at many phases of the pulsation, providing an unprecedented opportunity 
to test non-linear pulsation models. The full analysis of the PTI dataset will shed light 
on pulsation characteristics as a function of spectral type (O-rich, C-rich Miras) and 
wavelength, and we look forward to more of this work in the near future. 
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Figure 20. Pulsation curves for two different Mira Variables (at different 
wavelengths), (left panel) This panel shows the large amplitude pulsation of the S- 
Mira x Cyg seen at 905nm by COAST, although no obvious pulsations were seen at 
1290 nm (see Figure 2a in Young et al., 2000b, reproduced here with permission of 
Blackwell Publishing), (right panel) Here we see a well-sampled 2.2/xm pulsation curve 
(both open and filled plot symbols) for the O-rich Mira S Lac, from a recent PTI 
campaign (see Figure 2 in Thompson et al., 2002, reproduced with permission of the 
AAS). 



Optical Interferometry in Astronomy 



49 



4.1.4- Cepheid Pulsations Calibrate Period-Luminosity Relation Although phe- 
nomenologically related to measurements of pulsating AGB stars, observations of 
Cepheids are quite distinct in their scientific goals. As has been discussed in Sasselov 
and Karovska (1994) and earlier (e.g., Davis, 1976), optical interferometry will play an 
important role in independently calibrating the Cepheid distance scale, a crucial element 
of the cosmic distance ladder. By measuring the changing diameter of a nearby Cepheid 
and the coeval radial velocity curve through a pulsation cycle, the distance can be di- 
rectly inferred via the Baade-Wessehnk method. A flurry of initial results have appeared 
from GI2T (Mourard et al, 1997), NPOI (Armstrong et al, 2001), IOTA (Kervella et al., 
2001), and PTI (the first definitive detection of Cepheid pulsation; Lane et al., 2000b). 
However, most current published reports only weakly detect the pulsation, and defini- 
tive results will require longer baselines and/or much higher signal-to-noise ratio fringe 
measurements. 
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Figure 21. These data show the most recent pulsation curves for two Cepheid 

variables, rj Aql and Q Gem. The data was taken by the PTI interferometer at 1.65/im 
and, when combined with radial velocity data, result in the most accurate distances to 
these important primary distance indicators. This figure reproduced from Lane et al. 
(2002, see Figure 1), with permission from the AAS. 

This field is rapidly developing, both observationally and theoretically. PTI workers 
recently pubhshed a second Cepheid results paper (Lane et al., 2002), and these 
remarkable data are reprinted here in Figure 21. The derived Cepheid distances indeed 
agree with those measured from the Hipparcos parallax satellite, but the interferometer 
results are more precise. The distance uncertainty is now as (or more) affected by 
our uncertainty in physics of the stellar atmospheres; fortunately, Marengo et al. 
(2002) report timely theoretical studies of wavelength- and pulsational-phase-dependent 
variations in Cepheid limb-darkening and emphasize the need for careful calibration 
in order to interpret interferometry data accurately. With new observing campaigns 
underway at most interferometers, we can expect a rapid development here to really pin 
down the pulsational properties of Cepheids in the visible and near-infrared. In the near 
future, we can expect the best calibrated Cepheid distance scale to be interferometric. 



QGem 

1.9 I . . . 

1.85 
1.8 



„ 1.75 




0.75 0.25 0.5 0.75 1 0.25 
Phase 



Optical Interferometry in Astronomy 



50 



4.1.5. Imaging Stellar Surfaces As discussed in §3.1, the Cambridge group began 
interferometry research by using ( visible- hght) aperture masking on the WiUiam 
Hcrschel Telescope (WHT) in the Canary Islands (Baldwin et al., 1986b; Haniff et al., 
1987b) while developing the COAST interferometer. Interferometric imaging was 
performed and early results showed bright features (strong departure from circular 
symmetry) on the surface of Betelgeuse (Buscher et al., 1990), confirming some 
previous reports (e.g., Roddier and Roddier, 1983). No long-baseline (separate-element) 
interferometer would be able to investigate the nature of these features for years, and the 
Cambridge masking group has spent more than a decade since thoroughly investigating 
"hotspots" on red supcrgiants and giants. 

Over the last decade, it was shown that asymmetries are common (although 
not omnipresent) around red supergiants and giants at visible wavelengths (Wilson 
et al., 1992; Tuthill et al., 1997, 1999a), that these hotspots vary on a timescale of 
months (Wilson et al., 1997), and that the asymmetries become less-pronounced (even 
disappearing) into the infrared (Young et al., 2000a). The first image of a stellar 
photosphere using the COAST interferometer showed a featureless Betelgeuse (Burns 
et al., 1997). 

The hotspots were originally interpreted quite literally, as "hot" patches on 
the photosphere from upwellings of large convective elements (Schwarzschild, 1975). 
However, Young et al. (2000a) introduced a new paradigm which is illustrated in 
Figm^c 22. If the bright features were indeed caused by literal hotspots, one would not 
expect the features to complete disappear in the near-infrared. Here, we sec illustrated a 
model where the photosphere is surrounded by a molecular blanket (e.g., TiO) which is 
optically thick in the visible, but not in the infrared (^1/xm). Inhomogeneities (possibly 
caused by large-scale convection) allow visible hght to escape out of opacity holes. The 
results of Dyck and Nordgren (2002) indirectly support this model, by showing that 
even "continuum" visible diameters appear contaminated by TiO for late M-stars. 

4- 1.6. Binary Stars & Stellar Evolution Binary stars have been an indispensable tool 
for astronomers for centuries. Visual and spectroscopic observations yield reliable mass 
estimates and form the bedrock of stellar evolution theory (e.g., Eggen, 1967). The 
advent of speckle interferometry and optical long-baseline interferometry has led to a 
remarkable increase in the data quality and volume of binaries, including many short- 
period binaries for the first time (e.g., McAlister, 1985; Hartkopf et al., 2001). This work 
allows unprecedented tests of stellar evolution models on a case-by-case basis, through 
sub-1% precision of stellar parameters. 

In addition to the compelling science, the simple natiuc of their brightness 
distribution have made binary stars prime targets for most optical interferometers. 
Because of the high spatial resolution, these binaries tend to have short periods and 
thus full orbital elements can be determined by tracking the orbit. Not surprisingly, 
the first aperture synthesis images by the COAST and NPOI interferometers were of 
binary systems, and these first results are reproduced here in Figure 23. With sufficient 
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Figure 22. This figure shows a summary of the results from Young et al. (2000a), 
the most comprehensive investigation of the origin of surface hotspots on evolved 
stars. This synthesis shows nearly coeval images at different wavelengths of the 
surface of Bctclgcusc using a combination of WHT aperture masking and the COAST 
interferometer. Surface structures ("hotspots") apparent at visible wavelengths 
disappear when imaged in the infrared. The bottom-right panel shows the schematic 
model offered by these authors, where the hotspots are caused not by literally "hot" 
patches on the photospheric surface, but rather are caused by "opacity holes" in the 
molecular envelope (e.g., TiO) which allow visible light to escape in a patchy pattern. 
Figures appear here with permission of J. Young. 



angular resolution, interferometry yields the angular diameters of the components in 
addition to the binary separation vector and flux ratio. 

While earlier papers had concentrated on individual systems, Hummel et al. (1995) 
presented orbits of 8 systems with separations between 3 and 10 milliarcseconds using 
the Mark III interferometer. For some of the systems, precise mass and luminosity 
determinations allowed testing of stellar evolution models. 

In order to rigorously test stellar evolutionary models, the highest precision in 
parameters is needed. This requires combining data from multiple instruments and 
techniques in a global fit, and concomitant attention must be paid to systematic errors. 
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Figure 23. Baby Pictures: First true aperture synthesis images using long-baseline 
optical interferometry. (left panel) a. The binary star Capella reconstructed at 
two epochs using the COAST interferometer (see Figure 2 of Baldwin et al., 1996b, 
reproduced with permission of ESO). (right panel) Six epochs of Mizar A seen by the 
NPOI interferometer (see Figure 4 of Benson et al., 1997b, reproduced with permission 
of the AAS). 



Ideally, the orbital elements are fit directly to the visibility data and velocities, as 
implemented by Hummel et al. (1998) and Boden et al. (1999). In Figure 24, the 
results from a study of o Leo by Hummel et al. (2001) is presented. This study 
stands out because it combines interferometry data from the Mark HI, NPOI, and 
PTI interferometers, as well as radial velocity data, resulting in mass uncertainties of 
only ~0.5%. Stellar evolution isochrones can be put to a serious test for this system. 

While long-baseline interferometers allow very close binaries to be partially resolved 
(closest is probably 0.002" binary TZ Tri; Koresko et al., 1998) and wide systems to 
be characterized with incredible precision, this is not always very important. The most 
interesting science lies often in measuring unusual binary systems for the first time, 
such as the metal-poor double- lined binary system HD 195987 (Torres et al., 2002). 
Attractive targets for next generation of interferometer observations include systems 
with short-lived components such as Wolf-Rayet stars or Young Stellar Objects, since 
much less is already known about the masses of these objects. 

Lastly, the push for high dynamic range imaging of binary stars has obvious 
implications for detecting low-mass companions, even extrasolar planets, around nearby 
stars. This topic will be discussed further in §5.1.3 
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Figure 24. Precision binary parameters can be derived by combining interferometry 
and spectroscopy, as shown here for o Leo by Hummel et al. (2001). The left 
panel shows the derived orbit and data from multiple interferometers and the right 
panel shows a matching stellar isochrone along with the effective temperatures and 
luminosities of the two components (see Figures 10 & 11 in Hummel et al., 2001, 
reproduced with permission of the AAS). 



4-2. Circumstellar Environments 

Interferometers can also be used to probe the environments around stars, both at visible 
light and infrared wavelengths. As the sensitivity of facilities increase, lower surface 
brightness features can be measured, opening up new avenues of research. Advances in 
interferometric imaging are particularly relevant here, since gas and dust around stars 
might not be distributed uniformly and may be changing in time. While there has not 
been true imaging accomplished by long baseline interferometers in this area yet, we have 
included some of the unexpected recent results from Keck aperture masking (Tuthill 
et al., 2000c). While equally impressive results have also appeared using adaptive 
optics and speckle interferometry (in particular by the Weigelt group), we highlight 
the masking results because the observing methods and data reduction closely parallel 
that of optical interferometry and more truly reflect future capabilities; they directly 
motivate excitement in the potential of interferometric imaging with milli-arcsecond 
resolution and point in new scientific directions. 

4.2.1. Ha Envelopes around Hot Stars WhSie almost all the early visible 
interferometers focused on angular diameters and binary stars, an interesting exception 
was observations of the bright Ha line around Be stars. This emission was expected 
to be more extended, and thus more easily resolvable, than the tiny photosphere itself. 
The envelope of 7 Cas was first resolved by Thom et al. (1986) using the 12T, and 
Mourard et al. (1989) saw evidence for an envelope in rotation by inspecting multiple 
spectral channels across the line itself using the G12T. With a good range of baselines, 
the Mark 111 was able to detect definite asymmetries in 7 Cas and C, Tau (Quirrenbach 
et al., 1993a, 1994); for the latter, the maximum entropy method was used to visuahze 
the data as a "phase-less" image and this result is shown in Figure 25a. This data 
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lacked Fourier phase information, and thus can not quahfy as a true aperture synthesis 
image; nonetheless, the maximum entropy procedure provided an innovative and useful 
tool for visualizing this asymmetric envelope. 

The high spectral resolution of G12T later also uncovered asymmetric emission in 
these Be star envelopes in Ha (Stee et al., 1995), and observed in other lines too (Stee 
et al., 1998). Figure 25b shows results from Vakili et al. (1998) (see also Berio et al., 
1999) which proposed that the emission line region is very one-sided and time- variable. 
The asymmetry in this figure was derived from the behavior of the Fourier phases and 
amplitudes across the line profile; this amounts to a kind of phase-referencing where the 
continuum emission surrounding the Ha emission is used as a reference signal allowing 
the Fourier phases to be recovered. The origin of this "one-armed oscillation" could 
result from radiative effects, the presence of a companion, or dynamical behavior in a 
non-spherical potential; these interesting results should be confirmed and explored by 
the current generation of imaging interferometers. 




mjltiar^isecorjcts 

Figure 25. Spectral line observations in Ha around Be stars have found extended, 
asymmetric envelopes; the full potential of these kind of investigations remain 
untapped, (left panel) a. Here wc reproduce the "image" of the Ha envelope of 
Tau from Mark HI data (Figure 3 from Quirrenbach et al., 1994). (right panel) b. 
These schematic illustrations of the same envelope at two later dates were based on 
data from the GI2T interferometer (see Figure 4 from Vakili et al., 1998). Both figures 
are reproduced with permission of ESO. 



4-2.2. Accretion Disks and Young Stellar Objects There has been surprising and rapid 
progress in studies of how dense accretion disks evolve around pre-main sequence stars. 
Just five years ago, simple accretion scenarios incorporating passively-heated flared disks 
(e.g., Hillenbrand et al., 1992; Hartmann et al., 1993; Chiang and Goldreich, 1997) 
were widely accepted, adequate to explain the spectral energy distributions (SEDs) of 
most (low- mass) T Tauri stars and the higher- mass Herbig Ae/Be systems. However, 
recent observations with higher spatial resolution suggest a richer set of phenomena. 
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and has produced much excitement. Direct observational hnks are now even being 
made connecting the fields of star formation and planet formation, focusing on how 
accretion disks evolve into protoplanetary disks and finally to debris disks and planets. 
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Figure 26. (left panel) This figure shows the first Ilcrbig Ae/Be star visibility curve 
ever measured (see Figure 1 by Millan-Gabet et al., 1999b). The IOTA data of AB Aur 
at 1.6/xm (H-band) and 2.2/im (K'-band) indicated a much larger dust shell or disk 
than expected from accretion disk models (dash-dotted curves), (right panel) The first 
classical T Taiiri stars were observed by the PTI interferometer and also showed larger 
than expected infrared sizes. Here is the original data for T Tau itself, an important 
source but whose analysis is complicated by uncertain calibration due to the presence 
of a nearby companion in the interferometer field-of-view (sec Figure la by Akeson 
et al., 2000a). Both figures are reproduced with permission of the AAS. 



Infrared interferometry is playing an important role in elucidating the earliest stages 
of planetary formation by probing the density and temperature structure of the disks 
presumably before planets form. Malbet et al. (1998) reported the first resolved Young 
Stellar Object (YSO), using near- infrared data from the PTI interferometer; the source 
was FU Ori, a rare type of T Tauri whose emission is dominated by accretion luminosity, 
and the disk size was found to be roughly consistent with expectations. 

The first "normal" YSO to be resolved with an optical interferometer was the Herbig 
Ae/Be star AB Aur, and Millan-Gabet et al. (1999b) found the near-IR emission to be 
much larger than expected using the IOTA interferometer. These young massive stars 
have high luminosities, and thus were the brightest/easiest type of young star to study 
initially. These workers published a survey of 15 total Herbig Ae/Be stars which strongly 
reinforced the initial finding of large near-IR sizes (Millan-Gabet, 1999; Millan-Gabet 
et al., 2001). Also, they found no evidence for disk structures - intriguingly all the 
initial data were consistent with spherical distributions of dust. The much-awaited 
measurements of classical T Tauri stars would come soon thereafter from the Palomar 
Testbed Interferometer (Akeson ct al., 2000a, 2002). The few T Tauris measured also 
showed near-IR emission a few times larger than expected from accretion disk models. 
Figure 26 shows the visibility data from these important papers for AB Aur and T Tau. 

While infrared imaging of YSOs has only recently become possible with new 3+ 
telescope arrays (and nothing published yet), aperture masking on the Keck telescope 
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was able to resolve two of the brightest examples. Tuthill et al. (2001, 2002) imaged 
the bright source LkHa 101 and found a bright ring of emission, interpreting it as the 
hot dust at the inner edge of an accretion disk. While the size was larger than expected 
for a geometrically-thin, optically-thick disk model (e.g., Hillenbrand et al., 1992), the 
position of the inner edge was consistent with dust evaporating at temperatures above 
~1500 K when illuminated by direct stellar radiation (the standard paradigm for dust 
shells around evolved stars; see Rowan-Robinson and Harris, 1982; Dyck et al., 1984). 
Further this relation could explain the "large" sizes seen by IOTA and PTI. Figure 27a 
shows images of LkHa 101 and also of emission-line star MWC 349 by the Keck aperture 
maskers. We note that a subsequent LkHa 101 paper (Tuthill et al., 2002) included the 
first mid-IR measurements of a YSO disk size (using the ISI interferometer), and also 
reported unexpected changes in the LkHa 101 disk emission. 
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Figure 27. (left •panel) a. Aperture masking interferometry was used at the Keck 
telescope to create these 2.2/im images of disks around YSOs. The left panel shows the 
LkHa 101 dust disk with evidence of a central hole (Tuthill ct al., 2001), while the right 
panel reveals the MWC 349A disk viewed nearly edge-on (Danchi et al., 2001). (right 
•panel) b. Compendium of measured sizes of Herbig Ae/Be and T Tauri systems with 
comparison to theoretical dust sublimation radii (from Figure 1 of Monnier and Millan- 
Gabet, 2002, reproduced with permission of the AAS). Observed sizes are typically 
many times larger than expected from "classical" accretion disk models (dotted line). 
The number of YSO measurements will rapidly grow with the current size surveys at 
the Keck and VLT Interferometers. 



The first results of nulling interferometry observations of Herbig Ae/Be stars (see 
§3.5.7) have been recently published. By masking a single large aperture, Hinz et al. 
(2001b) found the spatial extent of three sources to be unresolved using a ~4m nulling 
baseline, an unexpected result. Interestingly, the same simple disk models which under- 
predicted the near-infrared sizes were shown to over-predict the mid-IR sizes. Clearly, 
more data is needed to determine what models are appropriate for disks around young 
stellar objects - none of the "standard" ones seem to work at near- or mid-infrared 
wavelengths when it comes to spatial observations (the spectral energy distributions 
can be fit by a number of models) . 
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Recently, Monnier and Millan-Gabet (2002) summarized the current literature of 
near-infrared disk sizes, and one of the figures appear here as Figure 27b. As indicated 
above, the most interesting fact is that nearly all the disks measured are many times 
larger than expected from "standard" disk models (geometrically-thin and optically- 
thick disks), indicating a large dust-free central cavity. A few theorists are incorporating 
these results into their models, and I recommend the papers of Natta et al. (2001) and 
DuUemond et al. (2001) for further discussion. Understanding how planetary systems 
eventually form out of these disks will not be possible until we are more certain of 
the initial physical conditions of material within a few AU of young stars. There is a 
survey already underway at the Keck Interferometer (and planned soon for the VLTI) 
to include many more T Tauri and Herbig Ae/Be stars, and we can look forward to 
more developments in this exciting area. 

4-2.3. Dust Shells & Molecules in Evolved Stars Long before the dusty disks around 
young stars could be observed, interferometry techniques were used to characterize dust 
shells around evolved stars. I briefly mentioned the early history of this work in the 
mid-infrared in §3.1, and the capabilities of 3 m telescopes were exploited in a series of 
near-infrared speckle measurements in the 1980s (see especially Dyck et al., 1984). 
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Figure 28. (left panel) This figure shows temporal changes in the IK Tau visibility 
curve probing the dust shell at ll.lS/im as shown by Hale et al. (1997, Figure 3). (right 
panel) These changes can be visualized through a maximum entropy reconstruction of 
the dust shell, reproduced here from Figure 4 of Hale et al. (1997). Both figures are 
reproduced with permission of the AAS. 



The most important paper on dust shells was published by Danchi et al. (1994), 
reporting mid-infrared observations of 13 evolved stars by the ISI interferometer. The 
data were fitted by radiative transfer models of dust shells, which strongly indicated 
that mass-loss appeared continuous around some stars and episodic around others. At 
the time, mass-loss on the AGB was still considered to be spherically-symmetric and 
continuous, and these results began to dislodge this simplified picture. 

Subsequent papers concentrated on detailed modeling of individual sources, as the 
evidence for non-uniform outfiows (Monnier et al., 1997) and deviations from spherical- 
symmetry (Lopez et al., 1997) accumulated. Hale et al. (1997) found time- variable 
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features in the visibility curve of the 0-rich mira IK Tau, and attributed it to moving 
dust shells expanding in the outflow. Figure 28 shows the visibility data and a maximum 
entropy reconstruction to help visualize the change in the dust shell morphology; as 
for the "image" of ( Tau by Quirrenbach et al. (1994), these radial profiles were 
reconstructed from phase-less data and so do not necessarily represent a true "image." 
The changes were consistent with the expected outflow speeds (derived from maser 
observations), and allowed an independent distance estimate to this source. 

The idea that mass-loss is not uniform in time nor necessarily spherically symmetric 
found conflrmation using other related techniques. At the same time, near-infrared 
speckle and aperture masking were being pursued on the new 8-m class telescopes. 
Imaging of the carbon star IRC +10216 revealed an incredibly inhomogeneous and 
asymmetric dust shell (Haniff and Buscher, 1998; Weigelt et al., 1998; Tuthill et al., 
2000c), and enough data has been collected to sec the dust structures evolve with time; 
I reproduce a figure from Tuthill et al. (2000c) which exploits the diffraction-limit of 
the world's largest telescope to image the details of this nebula (see Figure 29). For 
readers interested in near-infrared imaging of dust shells using speckle interferometry 
and aperture masking, I refer the reader to many recent papers by the Weigelt group 
and the Keck masking team, which will not be reviewed here. 

Another recent area of progress is the successful combination of high spectral 
resolution with high spatial resolution in the mid-infrared. A filtcrbank spectrometer 
was constructed and installed on the ISI interferometer allowing interferometry data to 
be collected on mid-infrared absorption lines (Monnicr et al., 2000c, a,b). Polyatomic 
molecules such as ammonia and silane form in dense outfiows of evolved stars, but 
at orders of magnitude greater abundance than expected (e.g., Keady and Ridgway, 
1993). This mystery was partially solved when the ISI found that the molecules form 
much further out in the flow than expected, possibly related (for the case of Silane) to 
the depletion of SiS onto grains (Bieging and Tafalla, 1993). This work highhghts the 
potential for studying cosmochemistry using interferometry. 

4-2. 4- Colliding Winds While no Wolf-Rayet dust shells have yet to be observed using 
a long-baseline interferometer, the discovery of pinwhecl nebulae around these stars by 
Tuthill et al. (1999b) and Monnier et al. (1999) using Keck aperture masking represents a 
new direction for the current generation of infrared imaging interferometers. An example 
of the spiral pattern seen in the visibihty data of WR 104 is shown in Figure 30; from 
this data (and the closure phases) images can be reconstructed showing a beautiful 
spiral structure which appears to rotate on the sky with a period ~243 days. These 
spinning, spiral dust shells result from dust formation at the interface of two colliding 
winds in a Wolf-Rayet and 0-star binary. 

These sources are usually obscured by local dust and thus are faint visible sources 
(V<;15). The faint V magnitudes make these targets difficult not only for adaptive optics 
systems, but also for all current interferometers (except the ISI) which use visible-light 
detectors for star tracking and tip-tilt correction. This limitation hopefully will be 



Optical Interferometry in Astronomy 



59 



aaci 
a» 

100 

S 100 
I ° 

A 

i -IDS 



^: 

li'l 1 ■! II. Iiii. 1 <l r'-. 


am m Ii'iiiiml 1 1"!: 




(ULih:98) 


I ^^^^^ 
L., I.^;...... .. ° 1 


p:^^^^^;' ""1"; 

= {Apr 93) ""'"'^.r 1: 



-3IX! 'ZOO "40) m ZOC KX) ^ -tM d 130 2<K> 3DII 

MlliaiQM4«ntt iriiiiuwcoridfr 
Oaitairs.(%olP9sl4vt S 3 4. ^ 7 10 15 30 70 



Figure 29. Shows temporal evolution of the inhomogeneous and clumpy dust shell 
around carbon star IRC +10216 as seen at 2.2/xm using Keck aperture masking (Figure 
1 from Tuthill et al., 2000b, reproduced with permission of the AAS). 



eliminated at some facilities by using an infrared star-tracker (for example, IOTA has 
plans along these lines). 

5. The Future 

5.1. Exciting Trends and Near-Future Science Potential 

As has been said before, now is a practical time for reviewing the major achievements of 
optical interferometry, for we are entering a new era boasting facilities with significantly 
greater sensitivity, angular resolution, spectral resolution, and wavelength-coverage. In 
this section, I will give my views of some of the new capabilities and the expected science 
returns. 

One important trend that must be bolstered is the inclusion of theorists and 
modellers in the observations and interpretations of interferometry data. In many 
areas, the interferometry observations are outstripping the ready tools for analysis. 
For example, the wavelength-dependent and time-dependent diameters of AGB stars 
require a combination of time-dependent hydrodynamical atmospheres and sophisticated 
radiative transfer codes, a problem very challenging even with today's supercomputers. 
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Figure 30. (left panels) a. These figures show the 2.2/xm visibihty curves, 1-D and 
2-D, of WR 104 observed by Keck aperture masking, (right panel) b. This figure 
shows four epochs of imaging of WR 104 showing a morphology consistent with an 
Archimedean spiral rotating at a uniform rate (from Monnier, 1999, Figure 13.6) 

Understanding the hotspots seen on the surfaces of stars will required 3-dimensional 
simulations of stellar convection. Accretion disk physics around young stars should 
include magnetic fields and demand thoughtful considerations of gas and dust physics 
in a 2-D or 3-D context. Dust production in coUiding winds is very poorly understood, 
and poses a formidable numerical simulation problem. While tackling these difficult 
physical problems will require the new high-resolution data from optical interferometers, 
it is also true that input from the modellers and theorists is needed to guide and suggest 
experiments and observing strategies. 

Another general comment is that increasing the angular resolution usually means 
probing ever decreasing physical scales. Since interferometers often probe scales smaller 
than an AU, significant changes in time are expected for even small characteristic 
velocities (~km/s). This poses both a risk and an opportunity: a risk since data 
must be taken rapidly and efficiently to accurately capture snapshots of ever-evolving 
and changing environs, and an opportunity to include dynamics and time-evolution 
into our models and understanding. Observing the dynamics of circumstellar and/or 
stellar environments allow new physics to be understood, physics that usually can 
not be unambiguously reconstructed from typical datascts. Thus, I hope that new 
dynamical information will break theoretical stalemates which paralyze a number of 
fields. Interferometers have the opportunity to revolutionize the way we think of the 
universe: from distant "frozen" images of the past, to a dynamic and engaging unfolding 
of the present. 
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5.1.1. New Long Baselines New long baselines will allow unprecedented high resolution 
measurements on select sources. With sub-milliarcsecond resolution, one can measure 
the diameters of "small" sources which have largely eluded current surveys, such as hot 
stars and nearby low-mass stars. Distortions in the photospheric shapes of rapidly 
rotating stars or binary stars in nearly Roche-lobe filling systems can be directly 
detected. Limb-darkening studies of important objects, such as Cepheids, can be 
accomplished to put the Cepheid distance scale of firm direct footing. Further, long 
baselines make a variety of exoplanet studies possible, such as directly detecting 51 
Peg B-hke planets ("Hot" jupiters) or resolving planetary transits across the stellar 
disk. The NPOI, CHARA, and SUSI interferometers will possess the longest baselines 
in the near-term, while future projects such as the MRO or 'GHANA (on Mauna Kca, 
HI) might someday extend the resolution below even 0.10 milliarcseconds with >1 km 
baselines. 

5.1.2. Imaging Imaging with optical interferometry is currently tedious at best, and 
can only investigate simple objects such as resolved photospheres or binary stars. The 6- 
telescope systems of NPOI and CHARA will soon possess the capability of (comparably) 
excellent "snapshot" coverage, allowing more complicated and higher dynamic range 
imaging of select targets. The CHARA array can not be reconfigured and hence will 
only image well targets with the appropriately-sized structures - for a maximum baseline 
of ~330 m at 1.65/im, the optimum size scale is a few milliarcseconds. The NPOI 
interferometer can be reconfigured to "fill-in" the (u,v)-plane completely over time, and 
be adjusted for individual sources to optimally measure the needed visibility and closure 
phase information. In the longer term, the auxiliary telescope array at the VLTI and 
the proposed outrigger telescopes at Keck will allow even fainter infrared targets to 
be observed. For imaging, the MRO is currently the most ambitious project in the 
works, hoping to include >10 telescopes, which would make it the premiere imaging 
interferometer in the world. 

Good imaging capabilities would open up new avenues of research, especially in 
studies of the circumstcllar environments at infrared wavelengths. The ability to study 
disks around young stars and the time evolution of gaps, rings, or other structures would 
revolutionize our understanding of planet formation. At visible wavelengths, imaging 
spots on the surfaces of other stars is a major goal, and would allow solar physics to be 
applied in detail to other stars for the first time. 

The unexpected discoveries of Keck aperture masking justify our optimism that 
imaging will uncover many new phenomena that currently are hidden unnoticed in 
spectral energy distributions. For example, the Wolf-Rayet dust spirals (see Figure 30) 
have only been observed in a few systems, and represent a new area of study when 
imaging intcrfcromctric arrays are fully commissioned. 

However, current imaging work using GOAST, NPOI, and IOTA interferometers 
suffer from the lack of dedicated software resources. Unfortunately, the decades of 
software development in radio interferometry can not be fully leveraged for optical 
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interferometry, since radio work now relies largely on phase-referencing techniques not 
generally available in the optical. New imaging software is needed which can take into 
account the unique nature of optical interferometry data as well as the different nature of 
our target sources. The recent adoption of a common data exchange format, defined by 
the COAST and NPOI interferometer teams, represent an important first step towards 
these goals (see http://www.mrao.cam.ac.uk/ jsylOOl/exchange/). 

5.1.3. Precision Interferometry This is a rapidly developing area since the advent 
of single-mode fibers for spatial filtering and "dual-star" phase referencing. When a 
model of the astronomical source is well-known, then incredibly precise measurements 
are possible. The most potential for this is in the general area of binary stars, where the 
stars either are point sources or partially-resolved uniform-disks. The case of detecting 
an exosolar planet is included in this category, since it can be considered as very high- 
dynamic-range imaging of a faint companion. 

While there are open questions in binary evolution and stellar astrophysics which 
demand such high precision, a more popular reason to pursue "Precision Interferometry" 
is towards detection of extrasolar planets around nearby stars. There are many ways 
this can be manifested, and I will outline a few of these. 

Narrow-angle astrometry is a comparatively "classical" way to detect an exosolar 
planet. Akin to the doppler shift-radial velocity method, precision astrometry attempts 
to detect the minute wobble of the parent star as a planet proceeds in its orbit. This can 
be done by monitoring the angular distance between a star and a background reference 
star. In this case, the target star is normally quite bright and used for phase-referencing 
to a faint star projected within an isoplanatic patch from the target (^30"). Lane et al. 
(2000a) reported the first measurements of this kind using the PTI (see Figure 31). 
For reference, the motion of Saturn and Jupiter perturb the Sun ~1 milliarcsecond as 
viewed from 10 pc. This technique will be applied by the Keck Interferometer and the 
VLTI interferometer for a planet survey, and there is talk of pursuing this in Antarctica 
where the isoplanatic patch is larger and the coherence times longer (e.g., Lloyd et al., 
2002; Swain, 2002). 

Another method also being aggressively pursued by the Keck and VLTI 
interferometers is a multi-wavelength approach to find massive exoplanets by detecting 
a very shght photocenter shift between different infrared bands due to hypothesized 
absorption bands in the planet's atmosphere (i.e., the differential phase method; e.g., 
Akeson and Swain, 1999; Lopez and Petrov, 2000). This method has the advantage 
of using the bright target star as its own phase reference. However, recent studies of 
line-of-sight variability of atmospheric water vapor (Akeson et al., 2000b) indicate that 
differential chromatic dispersion might be more difficult to calibrate for differential phase 
methods than originally expected. 

Precision measurements of closure phases can also be used to detect faint 
companions, a method which has not received as much attention. As described earlier in 
this review (§2.2.3), the closure phase is formed by summing the interferometer phases 
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Figure 31. State-of-the-art narrow-angle astrometry of the binary 61 Cyg by the PTI. 
For a period of one week, the residual astrometric error in declination was ~100 micro- 
arcseconds. Figure printed with permission of SPIE, originally appearing in Lane et al. 
(2000a). 



on three baselines around a triangle of telescopes, and this quantity is immune to 
atmospheric phase delays. The lack of attention to precision closure phase methods 
is understandable since few interferometers possess the requisite minimum of three 
telescopes. Monnier (2002) and Segransan (2002) recently discussed how closure phases 
are immune to dominant calibration problems of differential phase and that they can 
also be used to solve for all the parameters of a binary system without needing to 
measure any visibility amplitudes. For reference, a typical closure phase for a binary 
with brightness ratio of 10^ is ~0.01 degrees as long as the component separation is 
resolved by the interferometer - the same magnitude effect as for differential phase 
methods. 

Current published measurement precision of closure phases is only 0.5 to 5 degrees 
(Tuthill et al., 2000c; Benson et al., 1997b; Young et al., 2000a). Improving the three 
orders of magnitudes needed to detect even the brightest possible exoplanet is a daunting 
challenge. While there are surely unconsidered systematic effects (perhaps due to 
birefringence or drifts in optical alignment) which will degrade the sensitivity of the 
precision closure phase technique, the lack of any "showstopper" effects, like differential 
atmospheric dispersion for the differential phase methods, strongly argues for the further 
development of the closure phase technique. 

5.1.4- Nulling Another approach being pursued for planet detection is nulling 
(Bracewell, 1978). The initial nulling experiments with the MMT (Hinz et al., 1998) 
have continued (Hinz, 2001), and ultimately will be applied on the Large Binocular 
Telescope Interferometer (Hinz et al., 2001a). This project is still many years away, but 
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offer an alternative approacli to tlie "precision" pfiase metliods above. 

In tlie nearer term, tfie Keck Interferometer will be applying nulling in the mid- 
infrared in order to measure and characterize the zodiacal dust around nearby stars. 
This source of infrared radiation is expected to be the dominant background for an 
eventual space-based planet detection interferometer, the so-called "Terrestrial Planet 
Finder" mission (more information in §5.2). Serabyn and Colavita (2001) describe the 
"fully symmetric" nulling combiner being implemented on the Keck Interferometer, and 
initial on-sky tests are expected to begin in 2003. A more complete description of 
the observing strategy and expected sensitivity has been documented in Kuchner and 
Serabyn (2003). 

Nulling can also be applied on large single-apertures, and then are called nulling 
coronagraphs. New clever designs in coronagraphy are competing with nulling 
interferometry for space mission concepts to detect tcrrestial planets around other stars, 
and I recommend interesting papers on optimally-shaped and apodized pupils (Spergel, 
2002; Nisenson and Papaliohos, 2001), bandpass-limited image-masks (Kuchner and 
Traub, 2002), and phase-mask-based approaches (e.g., Guyon et al., 1999; Rouan et al., 
2000). 

5.1.5. Spectroscopy There have been only a few significant results combining 
spectroscopy and interferometry; fortunately, this is about to change. The near- 
infrared AMBER instrument, slated to arrive at the VLTI interferometer in 2003, will 
combine three telescope beams together and disperse the light with 3 different spectral 
resolutions, the maximum is R>10000. This resolution will allow interferometry on 
individual spectral lines in the 1-2.5/im regime, opening up shock-excited emission lines, 
CO-absorption/emission features, and even emission from YSO jets to be probed in 
novel and exciting ways for the first time. We can expect the value of interferometric 
observations to be greatly enhanced by these new capabilities. 

5.1.6. Polarimetry Imaging stars in polarized light with interferometers also promise 
fascinating new insights into many areas of astrophysics, although this capability 
is difficult to implement with current interferometers. Vakili et al. (2002) discuss 
interesting applications of combining the high spectral resolution of AMBER with 
polarimetry, and highlight the new capabilities for imaging scattered light and 
potentially even measuring stellar magnetic fields from the Zeeman effect. Experimental 
efforts (Rousselet-Perraut et al., 1997) in this area have been very hmited compared 
to the theoretical progress (Rousselet-Perraut et al., 2000); this situation should be 
remedied soon. 

5.1.7. New Observables Along with greater spectral coverage and more telescopes 
come new interferometric observables. While §5.1.3 discussed possible applications 
of differential phase and differential closure phase, there are other interferometric 
observables yet to be exploited for precision interferometry. 
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Measuring the diameter of a star by precisely locating the first null of the visibility 
pattern is immune to amplitude calibration errors. This could be done by using a 
well-calibrated spectrograph to search for the null, either measuring fringe amplitudes 
or looking for the signature phase-flip across the null (e.g., Mozurkewich, private 
communication). This technique is similar to the method of A. Michelson in measuring 
the diameter of Betelgeuse (Michelson and Pease, 1921), where the basehne was adjusted 
in order to find the visibility minimum as detected by his eyes. 

The closure amplitude (requires sets of 4 telescopes) is an important quantity in 
radio interferometry to compensate for unstable amplifier gains and varying antenna 
efficiencies that can be linked to individual telescopes (e.g., Readhead et al., 1980). 
Closure amplitudes are not practical for current optical interferometers partially because 
most fringe amplitude variations are not caused by telescope-specific gain changes 
but rather by changing coherence (e.g., due to changing atmosphere). However, the 
introduction of spatial filtering (e.g., single- mode fibers) should make the closure 
amphtude a useful tool for optical interferometry soon (see discussion in Monnier, 2000). 

Necessarily, most new observables have yet to be used in practice or described in 
print. I mention here a few possibilities that this author has considered to encourage 
future experimentation. For instance, it may be possible to use closure amplitudes in the 
case when fringe jitter causes loss of visibility contrast in a fringe-tracking interferometer, 
due to the way in which small random phase errors degrade coherence. Also, the Closure 
Differential Phase is a recently defined quantity (Monnier, 2002), introduced to overcome 
one limitation of current phase-referencing techniques, namely, that differential phase 
(and differential closure phase) methods requires assumptions about the source structure 
of the phase calibrator. 

5.1.8. Sensitivity (iKeck and VLTI) Another area where we expect immediate progress 
is in observing new classes of faint objects for the first time. The Keck and VLTI 

interferometers will have the capability of observing sources as faint as K~ll magnitude 
(down to K~20 with phase referencing), opening up extragalactic sources for the first 
time. By the time this review is printed, 1 expect that the first optical interferometric 
observations of the core of an Active Galactic Nuclei (AGN) will be announced. Size 
measurements of AGN should offer new constraints on models of the infrared continuum 
and, when coupled with high spectral resolution, could determine the physical origin of 
observed broad line regions and possibly even measure dynamical black hole masses. 

In terms of galactic sources, this increase in sensitivity will allow a broad census 
of sources to be taken, including YSOs spanning a broad range of ages, luminosities, 
and distances and binary systems of all masses. For instance, infrared observations of 
pre-main-sequence binaries allow unique probes of the evolution of binary fraction (e.g., 
Ghez et al., 1993) as well as important measurements of masses of young stellar objects 
(Tamazian et al., 2002). 1 expect interferometer observations to play an increasingly 
important role in this area as the sensitivity increases. 

Of course the additional sensitivity will permit new projects too, such as tracking 
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the motions of stars orbiting the black hole at the center of the Milky Way with an 
order of magnitude greater precision than possible today with single-apcrturc telescopes 
(e.g., Schodel et al., 2002). Precision astrometry may allow even new tests of General 
Relativity near super massive black holes at the center of nearby galaxies. 

In addition, the MIDI instrument for the VLTI will allow sensitive measurements in 
the mid-infrared for the first time. While the ISI interferometer pioneered interferometry 
in this wavelength range, MIDI+VLTI will be first to probe a wide range of sources 
with resolution of ~0."01 and down to N~4 mag (> 100 x fainter than the ISI). Mid- 
infrared observations are sensitive to emission from relatively cool dust and can peer 
through thicker layers of dust than possible in the visible or near-infrared. There are 
great possibilities for advancing our understanding of young and evolved stars both, and 
studying dust distributions in a variety of environments. 

5.2. Space Interferometry 

The greatest limitations to optical interferometers arise from atmospheric turbulence. It 
dramatically limits the sensitivity, the ability to do imaging, and forces the engineering 
to be clumsy and complicated. Space is naturally an ideal place for interferometry, with 
no atmosphere to corrupt the phase nor limit the coherent integration time. And long 
baselines are obviously possible by combining light intercepted by separate spacecraft 
fiying in formation. 

5.2.1. Critical Technologies Needed In order to successfully build space interferometers, 
many technologies must first be developed. To this day, there has not been any dedicated 
space interferometer fiown (except for the Fine Guidance Sensors on the Hubble Space 
Telescopes; e.g. Franz et al., 1991). 

For interferometers deployed on a single structure, one has to contend with truss 
vibrations, thermal and gravitational gradients, and an unusually large number of 
mechanisms (failures of which could end the mission). There are issues with propellant 
and power consumption for maneuvering the array to point around the sky. The Space 
Interferometry Mission (SIM) is in advanced planning stages and is being designed to 
measure accurate positions of stars with micro- arcsecond resolution. SIM is a "simple" 
2-element interferometer on a deployable truss (~10m maximum baseline), and will be 
the first space mission to attempt space interferometry. 

Ultimately, one would want to have baselines much longer than ~10 meters, and 
this will require separate, free-flying spacecraft. For a space interferometer consisting of 
"free-flyers," there are other problems. For instance, maintaining the physical distances 
between space telescopes to sub-micron tolerances is indeed a challenge. Probably this 
can not be done; however by monitoring the spacecraft drifts in real-time using laser 
metrology, the changing distances can be compensated for by onboard (short) delay lines. 
Some engineering missions have been proposed to test ideas, but have yet to really get- 
off-the-ground (e.g., the NASA Starhght mission was recently cancelled). NASA and 
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ESA should give such a test mission a high priority since the science potential for a 
free-flyer interferometer is so much greater than for one limited to a single structure. 

5.2.2. Review of Current NASA & ESA Missions There are a number of mission 
concepts involving space interferometry being considered by NASA and the European 
Space Agency (ESA). As mentioned before, the only one in advanced design stages 
is the NASA Space Interferometry Mission (SIM). In Table 4, I summarize some of 
the missions that are being proposed, and their main science drivers. Considering the 
unreliability of expected launch dates, I have omitted these from the table - it is unlikely 
any of these will fly before 2010 (2020?). 

NASA and ESA have spent much energy on designing missions to detect Earth- 
like planets around nearby stars, and to measure their crude reflectance (or emission) 
spectra. With luck, an extrasolar planet spectrum could encode distinctive atmospheric 
spectral features indicating the presence of life (biomarkers) on the distant planet (e.g., 
Woolf et al., 2002). While originally envisioned as an infrared interferometer mission, 
concepts involving a visible-light coronagraph have been proposed lately. This mission 
is known at the Terrestrial Planet Finder (TPF) in NASA, and as IRSI-Darwin at ESA. 
The summary table also includes a few TPF follow-on missions, such as "Life Finder." 
These missions are very futuristic, and testify to NASA's ebullient imagination. 

Another area of interest is imaging the far-infrared and sub-millimeter sky at high 
angular resolution using space interferometry. These wavelengths are difficult to access 
from the ground due to water absorption in the atmosphere. Because of this, the 
angular resolution of current observations are very limited (~30"); compared to all 
other wavelengths, the sky has been surveyed with the lowest resolution in the far-IR. 

The proposed NASA mission "Submillimeter Probe of the Evolution of Cosmic 
Structure" (SPECS) would be a separate-telescope space interferometer (possibly 
tethered together and not "free- flying" ) designed to map the sky with great sensitivity 
at a resolution comparable to that currently achievable at other wavelengths. (~0.010"). 
This would avoid the confusion-limited regime encountered by current low-angular- 
resolution galaxy count surveys, and allow the evolution of cosmic structure to be 
investigated back to high rcdshift. The SPIRIT mission is meant as a precursor to 
SPECS to test out various aspects on a single platform. 

The X-ray community has also proposed a space interferometer, which would boast 
micro-arcsecond resolution and be capable of studying the hot material at the event 
horizon of nearby Black Holes. Bolstered by successful lab experiments (Cash et al., 
2000), plans for a free-flying x-ray interferometer called the Micro- Arcsecond X-ray 
Imaging Mission (MAXIM) have begun. Controlling distances between macroscopic 
mirrors to picometer-precision, as is needed for X-ray interferometry, is indeed a 
daunting challenge. However, a MAXIM precursor mission with only a few meter 
baseline would have orders of magnitude greater resolution than the Chandra X-ray 
telescope and stands some chance of being flown. 
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Table 4. Proposed Space Interferometers 



Acronym 


Full Name & Primary Science Drivers 


QTTV/r /'AT A Q A ^ 

olM [SyAoA) 


Space Interferometry Mission 
Precision astrometry; exosolar planets 


T?1^0T /"NT A O A ^ 


Fourier-Kelvin Space Interferometer 

Find Jovian planets (nuUer); map circumstellar disks 


C!l\/rAX>T^ Q A ^ 

oMArti-o (HibAJ 


biVlAKi-o 

Test free-flying concept for ESA IRSI-Darwin mission 


TOOT T^«-«,,^i*^ ^T?C A A 

IHoi-Uarwm (JiioAJ 


Infra-Red Space Interferometer (one concept: Darwin) 
Image terrcstial planets (IR nuUer); measure spectra 


T'T^T? ^'\J A c A ^ 
i.rb (INAbAJ 


Terrestrial Planet Finder 

Image terrestial planets (IR nuUer); measure spectra 


CTDTOTT^ /AT A O A \ 
OFlKli (WADAj 


Space Infrared Interferometry Trailblazer 

Far-IR, sub-mm galaxy counts; precursor to SPECS 


CI3TrT^C? /AT \ Q W 

brJiOb (iNAbAj 


Submillimeter Probe of the Evolution of Cosmic Structure 
High-resolution map of High-Z universe (far-IR, sub-mm) 


SI (NASA) 


Stellar Imager 

Image surfaces of stars (visible, ultraviolet) 


MAXIM (NASA) 


Micro- Arcsccond X-ray Imaging Mission 

Alap l)lMck hole accretion disks and event horizons (X-rays) 


MAXIM Pathfinder (NASA) 


MAXIM Pathfinder 

Demonstrate feasibility of X-ray interferometry; achieve 100 /z-arcsecond resolution 


LF (NASA) 


Life Finder 

Search for biomarkers in planet spectra; TPF extension 


PI (NASA) 


Planet Imager 

Image surfaces of terrestial planets, 25x25 pixels 
(requires 6000km baselines, futuristic!) 



5.3. Future Ground-based Interferometers 

While it is interesting to speculate about the future of space interferometry, we recognize 
that it will be expensive, difficult, and slow-paced. In the next 10 or 20 years, we can 
expect more affordable and rapid progress to be possible from the ground. In this 
concluding section, I review some of the necessary characteristics of an Optical Very 
Large Array (OVLA). Ridgway (2000) discusses many of these considerations, and I 
refer the reader to his interesting report for further details. 

5.3.1. Design Goals The main design goal of a ncxt-gcncration optical interferometer 
array will be to allow the ordinary astronomer to observe a wide-range of targets without 
requiring extensive expert knowledge in interferometer observations. An imaging 
interferometer with great sensitivity could fulfill this promise by providing finished 
images, the most intuitive data format currently in use. It will not be a specialty 
instrument with narrow science drivers, but a general purpose facility to advance our 
understanding in a wide range of astrophysical areas. 
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5.3.2. Optical Very Large Array One way to achieve this design goal is to scale up the 
existing arrays. Simply put, this main goal will require an array with a large number 
of telescopes (>20 to allow reliable aperture synthesis imaging) and with large-aperture 
telescopes corrected by adaptive optics (preferably using laser guide stars for full-sky 
coverage), allowing a reasonably faint limiting magnitude (roughly speaking, brighter 
than ~15th magnitude in the infrared with no phase referencing). 

This array would hkely be reconfigurable, like the radio VLA, to allow different 
angular resolutions to be investigated. The longest baselines should cover a few 
kilometers (~0.1 milli-wccsecond resolution in the near-IR). The main limitation of such 
a system will be a small field-of-view, typically limited to the diffraction-limited beam 
of an individual telescope (for 10-m class telescopes, the instantaneous field of view 
would be only about ~50 milliarcseconds) - although mosaicing would be possible, as 
in the radio. There are schemes which can image a larger field simultaneously, but are 
probably not very practical. 

With an even larger (billion-dollar) budget, one can partially combine the goals of 
interferometry with the community priority for a 30 m diameter telescope. This clever 
idea was recently proposed by R. Angel and colleagues at the University of Arizona. 
In their "20/20" scheme, hght from two extremely large telescopes (diameter >20 
meters) would be combined in a Fizeau combination scheme, patterned after the Large 
Binocular Telescope, maintaining the entire field-of-view (~30", limited by atmospheric 
turbulence) with the rcsohition of the two-element interferometer. Further, this scheme 
maximizes raw collecting area and would boast potentially incredible sensitivity (>20 
mag!). One demanding feature of this design is that the two 20-|- m telescopes would 
have to smoothly move along a track in real-time to maintain the large field-of-view; 
this may not be impossible, but is surely an interesting complication. Further, the 
imaging advantages of this system only work when the two-telescope baseline is 5-10 x 
as large as the telescope diameter, and hence the "20/20" interferometer would have 
maximum baselines of only a few hundred meters at most, not much better than current 
interferometer arrays. While granting that this system could allow much fainter objects 
to be observed, this option would cost many times more than a dedicated OVLA system 
described above. 

5.3.3. Technological Obstacles Needed to be Overcome If optical interferometry is to 
continue its impressive growth over the coming decades, important breakthroughs must 
be made in critical areas. Here, I briefly list a few obvious improvements which would 
make an OVLA more affordable. 

The main advance needed to make the OVLA affordable will be the development 
of "cheap" large aperture telescopes with adaptive optics. Currently, it costs millions 
of dollars to build even a 4 m-class telescope - without adaptive optics. Advances in 
lightweight mirrors with adaptive optics designed-in from the beginning may change the 
economics of the situation. 

Another area which could revolutionize optical interferometry is advances in 
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photonic bandgap fiber materials (e.g., Mueller et al., 2002). These materials offer 
possibility of extremely wide-bandwidth, low dispersion and low-loss single-mode 
fibers, which could open up the possibility of practical fiber delay lines. Such an 
advance would greatly simplify the optical beam-train and engineering of an optical 
interferometer, making projects such as 'OHANA straightforward. This would put 
optical interferometry on more similar footing as radio interferometry, where cable delay 
lines (either coaxial or fiber) are routinely used. 

Combining dozens of telescopes may not be practical using bulk optics, and 
solutions involving integrated optics should be pursued. The main limitation of this 
technology is restricted wavelength coverage, currently only proven shortward of 2.2 /xm. 
Development of materials (e.g., lithium niobate) and fabrication processes that can 
extend the coverage into the thermal infrared (l-5;um) would mean that a general 
purpose interferometer could be built around an integrated optics combiner. Work is 
currently underway in Europe towards this end, in particular in pursuit of mid-infrared 
nulling capabihties for the ESA IRSI-Darwin mission (Haguenauer and others, 2002). 

Lastly, improved infrared detectors are crucial to maximizing the scientific output of 
a future interferometer. It has already been discussed here (see §3.3.5) that near-infrared 
detectors remain hmited by avoidable detector "read" noise, and a future OVLA must 
have better detectors. 

6. Conclusion 

After decades of development, optical interferometry is now poised to play a major role 
in mainstream astronomy. The emergence of well-funded interferometer "facilities," 
in particular the Very Large Telescope Interferometer and the Keck Interferometer, 
promise to revolutionize the impact of high-resolution observations in many areas of 
astrophysics. Clearly, the main beneficiaries will be stellar astrophysics and galactic 
astronomy, in particular the areas of star and planet formation, fundamental stellar 
properties, and all stages of stellar evolution. In addition, we can look forward to the 
first cxtragalactic results. 

Although the Keck and VLT interferometers hold immense promise, the field 
is currently driven forward by the activities of many other smaller groups, and 
scientific results will be dominated by these workers for the near future. While many 
experimental (astro)physics fields have matured to the point where future progress 
rests in "big science" collaborations and national research centers (e.g., NASA), optical 
interferometry represents one of the few healthy and active "experimental astrophysics" 
endeavors left in astronomy where university-based groups continue to make important 
technical innovations and astronomical discoveries. 

It is widely acknowledged that astronomy as a whole is experiencing a golden age 
of progress, spurred on by observational advances across the electromagnetic spectrum. 
Optical Interferometry has expanded in response to its own promising initial results, and 
we in the field look forward to exploiting the significant infrastructure buildup just now 
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being completed. I hope that the next review of optical interferometry will vindicate 
my optimism in the field, and that the pioneering discoveries reported here presage even 
grander exploits. It is safe to predict that the next decade will be critical to the field 
of high-resolution optical astronomy, since the scientific impact of current facilities will 
wholly determine whether the substantial funding required for an "Optical" Very Large 
Array can be justified to the international astronomical community. 
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